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INTRODUCTION 


The  release  of  certain  chemical  species  into  the  upper  atmosphere  results 
in  luminous  clouds  that  display  the  resonance  electronic-vibrational-rotation 
spectrum  of  the  released  species.  Such  spectra  are  seen  in  rocket  releases  of 
chemicals  for  upper  atmospheric  studies  and  upon  reentry  into  the  atmosphere  of 
artificial  satellites  and  missiles.  Of  particular  interest  in  this  connection 
is  the  observed  spectra  of  certain  metallic  oxides  and  air  diatomic  and  tri- 
atomic  species.  From  band  intensity  distribution  of  the  spectra,  and  knowledge 
of  the  f-values  for  electronic  and  vibrational  transitions,  the  local  conditions 
of  the  atmosphere  can  be  determined  (Ref.  1). 

Present  theoretical  efforts  which  are  directed  toward  a  more  complete  and 
realistic  analysis  of  the  transport  equations  governing  atmospheric  relaxation, 
including  chemical  effects,  and  the  propagation  of  artificial  disturbances  require 
detailed  information  on  atomic  and  ionic  reaction  rates  and  on  thermal  opacities 
and  LWIR  absorption  in  regions  of  temperature  and  pressure  where  molecular  effects 
are  important  (Refs.  2  and  i) .  Although  various  experimental  techniques  have  been 
employed  for  both  atomic  and  molecular  systems,  theoretical  studies  have  been 
largely  confined  to  an  analysis  of  the  properties  (bound-bound,  bound-free  and 
free-free)  of  atomic  systems  (Refs.  4  and  5).  This  has  been  due  in  large  part  to 
the  unavailability  of  reliable  wavef unctions  for  diatomic  molecular  systems,  and 
particularly  for  excited  states  or  states  of  open-shell  structures.  More  recently 
(Refs.  6-9),  reliable  theoretical  procedures  have  been  prescribed  for  such  systems 
which  have  resulted  in  the  development  of  practical  computational  programs. 

The  theoretical  analysis  of  atmospheric  reactions  requires  the  knowledge  of 
the  electronic  structure  of  atoms,  ions  and  small  molecular  clusters  of  nitrogen 
and  oxygen  and  the  interaction  of  water  or  other  small  molecules  with  these  clusters. 
Knowledge  of  the  chemistry  of  metal  oxide  species,  which  might  be  present  in  a 
contaminated  atmosphere,  is  also  desired.  In  addition  the  basic  collisional 
processes  involving  electrons,  ions  and  neutral  particles  must  be  understood  to 
evaluate  the  dynamic  effects  of  the  chemistry  of  the  atmosphere.  One  important 
reaction  is  the  charge  transfer  in  N+  +  0  ■>  N  +  0+.  This  has  been  studied 
experimentally  by  Neynaber  (Ref.  10  and  11)  at  collision  energies  between  0.5  eV 
and  25  eV  but  low  energy  (<  0.5  eV)  data  are  apparently  extremely  difficult  to 
measure . 

This  reaction  has  also  been  investigated  theoretically  by  Michels  (Ref.  12) 
who  identified  the  principle  states  (-*E+  and  ^n)  involved  in  this  collision  and 
has  carried  out  charge  transfer  calculations  of  the  cross-sections  for  this 
reaction  down  to  thermal  energies. 


As  an  extension  of  our  coliisional  studies  of  +  0,  an  examination  of 
Che  branching  ratios  in  the  dissociative-recombination  of  e  +  N0+  -+■  N(^D,^S)  +  0 
was  carried  out.  The  branching  into  N(^D)  and  N(^S)  has  been  shown  by  our  studies 
(Kef.  13)  to  be  governed  by  states  of  ^Z+  symmetry  of  NO.  Calculations  of  the 

n  a 

branching  from  an  initial  recombination  into  the  I  £  state  of  NO  to  atomic 
separations  of  N(^S)  4-  0(^P)  and  N(^D)  +  0  (^P)  have  been  carried  out.  This  study 
utilized  the  results  of  our  previous  research  on  the  electronic  structure  of 
N0+,  and  NO  (Technical  Reports  AFCRL-TR-75-0509  and  AFCRL-TR-76-0120) . 

Another  reaction  which  is  very  important  in  determining  the  properties  of 
the  ionosphere  is  the  ion-molecule  reaction: 

0+  +  N2  NO+  +  N 

This  reaction  provides  a  major  route  for  the  formation  of  N0+  and,  through 
dissociative-recombination  of  e  +  NO  ,  the  subsequent  production  of  N  in  the 
metastable  state.  Of  particular  interest  are  the  detailed  state-to-state  kinetic 
cross-sections  for  this  reaction  as  a  function  of  the  kinetic  energy  of  the  collision 
and  the  vibrational  state  of  No  in  the  ground  X  ''"l*  state.  The  rates  for  colli- 
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sions  with  metastable  N2  AJEU  are  also  of  importance  as  is  the  determination  of 
the  product  distribution  which  may  be  vibrationally  or  electronically  excited 
NO+,  or  N2+  +  0  at  sufficiently  high  collision  energies. 

Because  of  the  difficulty  of  conducting  experiments  to  measure  these  unknown 
cross-sections,  the  development  of  a  sound  theoretical  method  for  calculating  these 
state-to-state  processes  appears  necessary.  Although  relatively  little  work  of 
this  nature  has  been  done  in  the  past,  enough  theoretical  work  is  available  to 
indicate  that  the  development  of  such  procedures  can  be  made  practical,  particularly 
if  good  wave  functions  and  potential  energy  surfaces  are  available  for  the  inter¬ 
acting  species. 

The  actual  calculations  of  the  dynamics  of  molecular  collisions  are  ordinarily 
carried  out  with  the  aid  of  the  Born-Oppenheimer  separation  of  the  electronic  and 
nuclear  motion.  One  proceeds  by  calculating  the  electronic  energy  as  a  function  of 
the  positions  of  the  nuclei,  which  are  assumed  to  be  stationary.  This  electronic 
energy,  plus  the  electrostatic  repulsion  between  the  nuclei,  defines  a  potential 
energy  hypersurface  on  which  the  nuclei  may  be  regarded  as  moving.  A  potential 
energy  hypersurface  defined  in  this  way  is  referred  to  as  adiabatic,  and  is 
appropriate  for  describing  the  nuclear  motion  in  the  limit  of  low  velocity.  There 
are  many  collisions  for  which  an  adiabatic  potential  energy  hypersurface  provides 
an  adequate  description.  However,  most  reactive  collisions  and  many  collisions 
which  do  not  lead  to  reaction  are  inadequately  described  by  an  adiabatic  potential 
energy  hypersurface.  These  collisions  are  characterized  by  velocities  of  nuclear 
motion  sufficient  to  affect  adversely  the  Born-Oppenheimer  separation,  with  the 
result  that  the  overall  wavefunction  must  be  described  as  a  superposition  of  terms 
involving  different  electronic  energy  states.  Under  these  conditions,  it  will  be 
useful  to  consider  adiabatic  potent  i a  1  energy  hypersurfaces  corresponding  to  all 
electronic  states  relevant  to  the  overall  wavefunction. 
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When  the  different  potential  energy  hypersurfaces  are  well  separated  in 
energy,  the  nuclear  motion  can  ordinarily  be  described  in  terms  of  motion  on  a 
single  hypersurface.  However,  when  two  or  more  hypersurfaces  are  close  in  energy, 
they  can  be  expected  to  mix  appreciably  in  the  overall  wavefunction,  and  it  is 
then  necessary  not  only  to  calculate  the  hypersurfaces  but  also  to  calculate  the 
quantities  needed  to  discuss  their  mixing  in  the  overall  wavefunction. 

The  calculation  of  a  point  on  a  potential  energy  hypersurface  is  equivalent 
to  calculating  the  energy  of  a  diatomic  or  polyatomic  sysi_em  for  a  specified  nuclear 
configuration,  and  therefore  will  present  considerable  practical  computational 
difficulty.  For  certain  problems  or  nuclear  configurations,  the  maximum  possible 
accuracy  is  needed  and  under  these  conditions  relatively  elaborate  ab  initio 
methods  are  indicated.  Under  other  conditions,  it  may  be  possible  to  use  less 
elaborate  and  more  rapid  computational  methods,  and  density  functional  or  other 
approaches  may  then  prove  useful. 

The  research  program  contained  herein  was  devoted  to  a  theoretical  study  of 
the  energetics  and  kinetic  reaction  rates  of  0+  +  N2  collisions.  The  goal  of 
this  research  program  was  to  develop  technical  information  concerning  this  system 
relevant  to  DNA  interests  in  upper  atmospheric  reactions. 
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CURRENT  STATUS  OF  QUANTUM  MECHANICAL  METHODS  FOR 
DIATOMIC  AND  POLYATOMIC  SYSTEMS 


The  application  of  quantum  mechanical  methods  to  the  pr?diction  of 
electronic  structure  has  yielded  much  detailed  information  about  atomic  and 
molecular  properties  (Ref.  7).  Particularly  in  the  past  few  years,  the  avail¬ 
ability  of  high-speed  computers  with  large  storage  capacities  has  made  it 
possible  to  examine  both  atomic  and  molecular  systems  using  an  ab  initio 
approach,  wherein  no  empirical  parameters  are  employed  (Ref.  14).  Ab  initio 
calculations  for  diatomic  molecules  employ  a  Hamiltonian  based  on  the  non- 
relativistic  electrostatic  interaction  of  the  nuclei  and  electrons,  and  a 
wavefunction  formed  by  antisymmetrizing  a  suitable  many-electron  function  of 
spatial  and  spin  coordinates.  For  most  applications  it  is  also  necessary  that 
the  wavefunction  represent  a  particular  spin  eigenstate  and  that  it  have 
appropriate  geometrical  symmetry.  Nearly  all  the  calculations  performed  to 
date  are  based  on  the  use  of  one-electron  orbitals  and  are  of  two  types: 
Hartree-Fock  or  configuration  interaction  (Ref.  8). 

Hartree-Fock  calculations  are  based  on  a  single  assignment  of  electrons 
to  spatial  orbitals,  following  which  the  spatial  orbitals  are  optimized, 
usually  subject  to  certain  restrictions.  Almost  all  Hartree-Fock  calculations 
have  been  subject  to  the  assumption  that  the  diatomic  spatial  orbitals  are 
all  doubly  occupied,  as  nearly  as  possible,  and  are  all  of  definite  geometri¬ 
cal  symmetry.  These  restrictions  define  the  conventional,  or  restricted, 
Hartree-Fock  (RHF)  method  (Refs.  15  and  16).  RHF  calculations  can  be  made 
with  relatively  large  Slater-type  orbital  (STO)  basis  sets  for  diatomic  mole¬ 
cules  with  first  or  second-row  atoms,  and  the  results  are  convergent  in  the 
sense  that  they  are  insensitive  to  basis  enlargement.  The  RHF  model  is  ade¬ 
quate  to  give  a  qualitatively  correct  description  of  the  electron  interaction 
in  many  systems,  and  in  favorable  cases  can  yield  equilibrium  interatomic 
separations  and  force  constants.  However,  the  double-occupancy  restriction 
makes  the  RHF  method  inappropriate  in  a  number  of  circumstances  of  practical 
interest.  In  particular,  it  cannot  provide  potential  curves  for  molecules 
dissociating  into  odd-electron  atoms  (e.g.,  NO  at  large  internuclear  separa¬ 
tion),  or  into  atoms  having  less  electron  pairing  than  the  original  molecule 
[e.g.,  O2  ■*  0 ( ^ P) ] ;  it  cannot  handle  excited  states  having  unpaired 
electrons  (e.g.,  the  ^ T.  states  of  O2  responsible  for  the  Schumann-Runge 
bands);  and,  in  general,  it  gives  misleading  results  for  molecules  in  which 
the  extent  of  electron  correlation  changes  with  internuclear  separation. 

Configuration-interaction  (Cl)  methods  have  the  capability  of  avoiding 
the  limitations  of  the  RHF  calculations.  If  configurations  not  restricted  to 
doubly-occupied  orbitals  are  included,  a  Cl  can,  in  principle,  converge  to 
an  exact  wave-function  for  the  customary  Hamiltonian.  However,  many  Cl 


calculations  have  in  fact  been  based  on  a  restriction  to  doubly-occupied 
orbitals  and  therefore  retain  many  of  the  disadvantages  of  the  RHK  method 
(Ref.  8).  The  use  of  general  Cl  formulations  involves  three  considerations, 
all  of  which  have  been  satisfactorily  investigated:  the  choice  of  basis 
orbitals,  the  choice  of  configurations  (sets  of  orbital  assignments),  and 
the  specific  calculations  needed  to  make  wavefunct ions  describing  pure  spin 
states  (Ref.  6).  The  first  consideration  is  the  art  associated  with  quantum 
mechanical  electronic  structure  calculations.  Many  mi  mods  (iterative  NSO, 
perturbation  selection,  first  order  Cl,  etc.)  have  been  advocated  for  the 
optimum  choice  of  configurations.  There  are  not  firm  rules  at  present  and 
the  optimum  choice  is  a  strong  function  of  the  insight  of  the  particular 
research  investigator.  The  last  consideration,  proper  spin  and  symmetry  pro¬ 
jection,  has  proved  difficult  to  implement,  but  computer  programs  have  been 
developed  for  linear  projection  algebra  at  this  Center,  and  the  Cl  method 
has  been  found  of  demonstrable  value  in  handling  excited  states  and  dissoc¬ 
iation  processes  which  cannot  be  treated  with.  RHF  techniques. 

Either  of  the  above  described  methods  for  ab  initio  calculations  reduces 
in  practice  to  a  series  of  steps,  the  most  important  of  which  are  the  evalua¬ 
tion  of  molecular  integrals,  the  construction  of  matrix  elements  of  the 
Hamiltonian,  and  the  optimization  of  molecular  orbitals  (RHF)  or  configuration 
coefficients  (Cl).  For  diatomic  molecules,  these  steps  are  all  comparable 
in  their  computing  time,  so  that  a  point  has  been  reached  where  there  is  no 
longer  any  one  bottleneck  determining  computation  speed.  In  short,  the 
integral  evaluation  involves  the  use  of  ellipsoidal  coordinates  and  the  intro¬ 
duction  of  the  Neumann  expansion  for  the  interelectronic  repulsion  potential 
(Ref.  17);  the  matrix  element  construction  depends  upon  an  analysis  of  the 
algebra  of  spin  eigenfunctions  (Ref.  18),  and  the  orbital  or  configuration 
optimization  can  be  carried  out  by  eigenvalue  techniques  (Refs.  19,  20).  All 
the  steps  have  by  now  become  relatively  standard  and  can  be  performed 
efficiently  on  a  computer  having  65,000  to  130,000  words  of  core  storage,  a 
cycle  time  in  the  microsecond  range,  and  several  hundred  thousand  words  of 
peripheral  storage. 

Both  the  RHF  and  Cl  methods  yield  electronic  wavefunctions  and  energies 
as  a  function  of  the  internuclear  separation,  the  RHF  methods  for  one  state, 
and  the  Cl  method  for  all  states  considered.  The  electronic  energies  can  be 
regarded  as  potential  curves,  from  which  may  be  deduced  equilibrium  inter¬ 
nuclear  separations,  dissociation  energies,  and  constants  describing  vibra¬ 
tional  and  rotational  motion  (including  anharmonic  and  rotation-vibration 
effects).  It  is  also  possible  to  solve  the  Schrodinger  equation  for  the  mo¬ 
tion  of  the  nuclei  subject  to  the  potential  curves,  to  obtain  vibrational 
wavefunctions  for  use  in  transition  probability  calculations.  The  electronic 
wavefunctions  themselves  can  be  used  to  estimate  dipole  moments  of  individual 
electronic  states,  transition  moments  between  different  electronic  states, 
and  other  properties.  While  all  of  the  calculations  described  In  this 


paragraph  have  been  carried  out  on  some  systems,  the  unavailability  of  good 

tlTt  TJh  WaVe£unCtions  and  Potential  curves  has  limited  actual  studies  of 
most  of  these  properties  to  a  very  small  number  of  molecules. 
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METHOD  OF  APPROACH 


Central  to  these  studies  are  the  actual  quantum-mechanical  calculations 
which  must  be  carried  out  for  the  molecular  species.  For  added  clarity, 
various  aspects  of  these  calculations  are  discussed  in  individual  subsections. 


1.  Levels  of  Approximation 

Much  evidence  on  diatomic  and  polyatomic  systems  indicates  the  near  ade¬ 
quacy  of  a  minimum  Slater-type-orbital  (STO)  basis  for  constructing  qualita¬ 
tively  correct  molecular  wavefunctions  (Refs.  21  and  22).  This  means  inner- 
shell  and  valence-shell  STO's  of  quantum  numbers  appropriate  to  the  atoms 
(Is,  2s,  2p  for  C,  N,  0;  Is  for  H) .  The  main  deficiency  of  the  minimum  basis 
set  is  in  its  inability  to  properly  describe  polarization  and  the  change  of 
orbital  shape  for  systems  which  exhibit  large  charge  transfer  effects.  Values 
of  the  screening  parameters  r,  for  each  orbital  can  either  be  set  from  atomic 
studies  or  optimized  in  the  molecule;  the  latter  approach  is  indicated  for 
studies  of  higher  precision.  When  high  chemical  accuracy  is  required,  as  for 
the  detailed  studies  of  the  ground  or  a  particular  excited  state  of  a  system, 
a  more  extended  basis  must  be  used.  Double-zeta  plus  polarization  functions 
or  optimized  MO's  are  required  for  reliable  calculated  results  of  chemical 
accuracy. 

The  chosen  basis  sets  give  good  results  only  when  used  in  a  maximally 
flexible  manner.  This  implies  the  construction  of  Cl  wavefunctions  with  all 
kinds  of  possible  orbital  occupancies,  so  that  the  correlation  of  electrons 
into  overall  states  can  adjust  to  an  optimum  form  at  each  geometrical  confor¬ 
mation  and  for  each  state.  Except  when  well-defined  pairings  exist,  as  for 
closed  shell  and  exchange  dominated  systems,  a  single-configuration  study 
(even  of  Hartree-Fock  quality)  will  be  inadequate. 


2.  Spin  and  Symmetry 

Proper  electronic  states  for  systems  composed  of  light  atoms  should 

2 

possess  definite  eigenvalues  of  the  spin  operator  S  as  well  as  an  appropriate 
geometrical  symmetry.  The  geometrical  symmetry  can  be  controlled  by  the 
assignment  of  orbitals  to  each  configuration,  but  the  spin  state  much  be  ob¬ 
tained  by  a  constructive  or  projective  technique.  Formulas  have  been  developed 
(Ref.  18)  for  projective  construction  of  spin  states  from  orthogonal  orbitals, 
and  programs  implementing  these  formulas  have  for  several  years  been  in  rou¬ 
tine  use  at  UTRC. 
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One  of  the  least  widely  appreciated  aspects  of  the  spin-projection 
problem  is  that  the  same  set  of  occupied  spatial  orbitals  can  sometimes  be 
coupled  to  give  more  than  one  overall  state  of  given  S  quantum  number.  It  is 
necessary  to  include  in  calculations  all  such  spin  couplings,  as  the  optimum 
coupling  will  continuously  change  with  changes  in  the  molecular  conformation. 
This  is  especially  important  in  describing  degenerate  or  near-degenerate 
excited  electronic  states. 


3.  Method  of  Ab  Initio  Calculation 

A  spin-free,  nonrelativistic ,  electrostatic  Hamiltonian  is  employed  in 
the  Born-Oppenheimer  approximation.  In  systems  containing  atoms  as  heavy  as 
N  or  0  this  approximation  is  quite  good  for  low-lying  molecular  states.  For 
a  diatomic  molecule  containing  n  electrons,  the  approximation  leads  to  an 
electrostatic  Hamiltonian  depending  parametrically  on  the  internuclear  separa¬ 
tion,  R: 


J(R) 


+ 


Za  zb 
R 


+ 


I 


i>j>l 


(1) 


where  ZA  and  Zg  are  the  charges  of  nuclei  A  and  B,  and  riA  is  the  separation 
of  electron  i  and  nucleus  A.  Jf  is  in  atomic  units  (energy  in  Hartrees, 
length  in  Bohrs) . 

Electronic  wavefunct ions  iJ/(R)  are  made  to  be  optimum  approximations  to 
solutions,  for  a  given  R,  of  the  Schrodinger  equation 


(2) 


by  invoking  the  variational  principle 


SW  (R)=  S 


fy'(R)*(R)c(T 


(3) 


The  integrations  in  Eq.  (3)  are  over  all  electronic  coordinates  and  the  sta¬ 
tionary  values  of  W(R)  are  approximations  to  the  energies  of  states  described 
by  the  corresponding  v|>(R) .  States  of  a  particular  symmetry  are  studied  by 
restricting  the  electronic  wavefunction  to  be  a  projection  of  the  appropriate 
angular  momentum  and  spin  operators.  Excited  electronic  states  corresponding 
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to  a  particular  symmetry  are  handled  by  construction  of  configuration- 
interaction  wavefunct ions  of  appropriate  size  and  form. 

The  specific  form  for  i*'(R)  may  be  written 


'Mr )  =  £  c (r) 


(4) 


where  each  liiy(R)  is  referred  to  as  a  configuration,  and  lias  the  general 
structure 


V(R )--A0s  n  4>  ,  (Li,R)0M 


1=1 


(r>) 


where  each  <j>^  is  a  spatial  orbital,  H  is  the  antisymmetrizing  operator, 

0S  is  the  spin-projection  operator  for  spin  quantum  number  S,  and  is  a 
product  of  a  and  8  ot.e-electron  spin  functions  of  magnetic  quantum  number  M. 

No  requirement  is  imposed  as  to  the  double  occupancy  of  the  spatial  orbital, 
so  Eqs.  (4)  and  (5)  can  describe  a  completely  general  wavefunct ion. 

In  Hartree-Fock  calculations  iJi(R)  is  restricted  to  a  single  which  is 
assumed  to  consist  as  nearly  as  possible  of  doubly-occupied  orbitals.  The 
orbitals  4^  are  then  selected  to  be  the  linear  combinations  of  basis  orbitals 
best  satisfying  Eq.  (3).  Writing 


~  Z  °j/i  X  v 


(6) 


the  av^  are  determined  by  solving  the  matrix  Hartree-Fock  equations 


l  FX„a„i  =  «iZ  Sx„a„j  (each  X)  (7) 

where  is  the  orbital  energy  of  4>u^. 

The  Fock  operator  F^^  has  been  thoroughly  discussed  in  the  literature 
(Ref.  16)  and  depends  upon  one-  and  two-electron  molecular  integrals  and  upon 
the  aV£.  This  makes  Eq.  (7)  nonlinear  and  it  is  therefore  solved  iteratively. 
UTRC  has  developed  programs  for  solving  Eq.  (7)  for  both  closed  and  open-shell 
systems,  using  basis  sets  consisting  of  Slater-type  atomic  orbitals.  Examples 
of  their  use  are  in  the  literature  (Ref.  6). 


q 


In  configuration  interaction  calculations,  the  summation  in  Eq.  (4)  has 
more  than  one  term,  and  the  £  are  determined  by  imposing  Eq.  (3),  to  obtain 
the  secular  equation 


X  ( — 


w 


AH' )  CM  =  0 


(each  p  ) 


(8) 


where 


h  aW*;  (R)J  (WvWc/t 

S/ii/  =f  *1  {R)*„  (r)  dT 


(9) 


Equation  (8)  is  solved  by  matrix  diagonalization  using  either  a  modified 
Givens  method  (Ref.  19)  or  a  method  due  to  Shavitt  (Ref.  20). 

The  matrix  elements  Hys;  and  may  be  reduced  by  appropriate  operator 
algebra  to  the  forms 


H/xy'Z  «P<0M  P 


,  R ) 


J  (R)P 


n 

n 

i  =  l 


(I 


'>> 


(10) 


0zp 


n  N 

n  (Xi > ^  \ 


(id 


where  P  is  a  permutation  and  Ep  its  parity.  The  sum  is  over  all  permutations. 
<0Mj0sP|0M>  is  a  "Sanibel  coefficient"  and  the  remaining  factors  are  spatial 
integrals  which  can  be  factored  into  one-  and  two-electron  integrals.  If  the 
are  orthonormal,  Eqs.  (10)  and  (11)  become  more  tractable  and  the  H  and 
S  may  be  evaluated  by  explicit  methods  given  in  the  literature  (Ref.  18). 
Computer  programs  have  been  developed  for  carrying  out  this  procedure,  and 
they  have  been  used  for  problems  containing  up  to  64  total  electrons,  10  un¬ 
paired  electrons,  and  several  thousand  configurations. 

The  Cl  studies  described  above  can  be  carried  out  for  any  orthonormal 
set  of  |  for  which  the  molecular  integrals  can  be  calculated.  Programs 
developed  by  UTRC  make  specific  provision  for  the  choice  of  the  $  ^  as 
Slater-type  atomic  orbitals,  as  symmetry  molecular  orbitals,  as  Hartree-Fock 
orbitals,  or  as  more  arbitrary  combinations  of  atomic  orbitals. 
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4.  Molecular  Integrals 


The  one-  and  two-electron  Integrals  needed  for  the  above  described  method 
of  calculation  are  evaluated  for  STO’s  by  methods  developed  by  the  present 
Investigators  (Ref.  23).  All  needed  computer  programs  have  been  developed 
and  fully  tested  at  UTRC. 


5.  Configuration  Selection 

Using  a  minimum  basis  plus  polarization  set  of  one-electron  functions,  a 
typical  system  can  have  of  the  order  of  10^  configurations  in  full  Cl  (that 
resulting  from  all  possible  orbital  occupancies) .  It  is  therefore  essential 
to  identify  and  use  the  configurations  describing  the  significant  part  of  the 
wavefunction.  There  are  several  way;;  to  accomplish  this  objective.  First, 
one  may  screen  atomic-orbital  occupancies  to  eliminate  those  with  excessive 
formal  charge.  Alternatively,  in  a  molecular-orbital  framework  one  may 
eliminate  configurations  with  excessive  numbers  of  anti-bonding  orbitals. 

A  third  possibility  is  to  carry  out  an  initial  screening  of  configurations, 
rejecting  those  whose  diagonal  energies  and  interaction  maLrix  elements  do  not 
satisfy  significance  criteria.  Programs  to  sort  configurations  on  all  the 
above  criteria  are  available  at  UTRC. 

Other  potentially  more  elegant  methods  of  configuration  choice  involve 
formal  approaches  based  on  natural-orbital  (Ref.  24)  or  multiconfiguration 
SCF  (Ref.  25)  concepts.  To  implement  the  natural-orbital  approach,  an  initial 
limited-CI  wave-function  is  transformed  to  natural-orbital  form,  and  the  re¬ 
sulting  natural  orbitals  are  used  to  form  a  new  Cl.  The  hoped-for  result  is 
a  concentration  of  the  bulk  of  the  Cl  wavefunction  into  a  smaller  number  of 
significant  terms.  The  multiconfiguration  SCF  approach  is  more  cumbersome, 
but  in  principle  more  effective.  It  yields  the  optimum  orbital  choice  for 
a  preselected  set  of  configurations.  This  approach  works  well  when  a  small 
number  of  dominant  configurations  can  be  readily  identified. 

It  should  be  emphasized  that  the  problem  of  configuration  choice  is  not 
trivial,  and  represents  an  area  of  detailed  study  in  this  research.  The 
existence  of  this  problem  causes  integral  evaluation  to  be  far  from  a  unique 
limiting  factor  in  the  work. 
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6.  Density  Function  Approach  -  Xa  Model 

The  Xa  model  (Ref.  26)  for  the  electronic  structure  of  atoms,  molecules, 
clusters  and  solids  is  a  local  potential  model  obtained  by  making  a  simple 
approximation  to  the  exchange  -  correlation  energy.  If  we  assume  a  non- 
relativistic  Hamiltonian  with  only  electrostatic  interactions,  it  can  be  shown 
that  the  total  energy  E  of  a  system  can  be  written  exactly  (Ref.  27)  (in 
atomic  units)  as 


E=  I  Oj  <  u,  -4“  +  I  T~~  ui  >  +  £  ~t 

.  i.  rl  il  I  li  1/  ' 


Vf  ¥ 


fj.  'IH- 


(12) 


This  expression  is  exact  provided  the  u^  are  natural  orbitals  and  n^  are  their 
occupation  numbers  (i.e.,  eigenfunctions  and  eigenvalues  of  the  first  order 
density  matrix).  The  first  term  in  Eq.  (12)  represents  the  kinetic  and  elec¬ 
tron-nuclear  energies.  The  second  term  is  the  nuclear  repulsion  energy.  The 
sums  (u,v)  are  over  all  the  nuclear  charges  in  the  system.  The  third  term  is 
the  electron-electron  repulsion  term,  which  represents  the  classical  electro¬ 
static  energy  of  the  charge  density  p  interacting  with  itself,  where 


*>(i)  *1  n.u*(D  u.(D  (13) 

i  1 

The  last  term  Exc  represents  the  exchange  correlation  energy  and  can  be  ex¬ 
pressed  formally  as 


E,c=  t/pU)  dr!  I  drT  ,  (14) 

12 

where  pxc  (1,2)  represents  the  exchange-correlation  hole  around  an  electron 
at  position  1.  In  the  exact  expression,  pxc  is  dependent  on  the  second-order 
density  matrix.  In  the  Hartree-Fock  approximation  Exc  is  the  exchange  energy, 
Pxc  represents  the  Fermi  hole  due  to  the  exclusion  principle  and  depends  only 
on  the  first  order  density  matrix.  In  the  Xa  method,  we  make  a  simpler 
assumption  about  pxc.  If  we  assume  that  the  exchange-correlation  hole  is 
centered  on  the  electron  and  is  spherically  symmetric,  it  can  be  shown  that 
the  exchange-correlation  potential 
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(15) 


is  inversely  proportional  to  the  range  of  the  hole,  rs,  where  r  is  defined 
by 


— T“  '1  =  l  (16) 

Therefore,  in  the  Xa  model,  the  potential  Uxc  is  proportional  to  p1/J(r). 

We  define  a  scaling  parameter  a  such  that 


U*<y°,=  ‘  ■¥"  (3/>0)/8ir)'/3  (17) 


The  expression  in  Eq.  (17)  is  defined  so  that  a  =  2/3  for  the  case  of  a 
free  electron  gas  in  the  Hartree-Fock  model  (Ref.  28)  and  a  =  1  for  the  po¬ 
tential  originally  suggested  by  Slater  (Ref.  29).  A  convenient  way  to  choose 
this  parameter  for  molecular  and  solid  state  applications  is  to  optimize 
the  solutions  to  the  Xa  equations  in  the  atomic  limit.  Schwarz  (Ref.  30)  has 
done  this  for  atoms  from  Z  x  1  to  Z  =  41  and  found  values  between  2/3  and  1. 

In  the  "spin  polarized"  version  of  the  Xa  theory,  it  is  assumed  (as  in 
the  spin-unrestricted  Hartree-Fock  model)  that  electrons  interact  only  with 
a  potential  determined  by  the  charge  density  of  the  same  spin.  In  this  case 
the  contribution  to  the  total  energy  is  summed  over  the  two  spins,  s  =  -  1/2. 

E*c  1  T  £  <l8) 

where  the  potential  is  spin-dependent 

U.  (3o(|)/4tt),/3  (19) 

*  O' « s  *  hs 


and  pg  is  the  charge  density  corresponding  the  electrons  of  spin  s.  The  spin- 
polarized  Xa  model  is  useful  for  describing  atoms  and  molecules  with  open- 
shell  configurations  and  crystals  which  are  ferromagnetic  or  anti-ferromagnetic. 
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Once  one  has  made  the  Xa  approximation  to  the  total  energy  functional  E 
in  Eq.  (12),  then  the  rest  of  the  theory  follows  from  the  application  of 
variational  principle.  The  orbitals  uj_  are  determined  by  demanding  that  E  be 
stationary  with  respect  to  variations  in  u^.  This  leads  to  the  set  of  one- 
electron  Xa  equations 


where  is  the  one-electron  eigenvalue  associated  with  u^.  Since  p(r)  is 
defined  in  terms  of  the  orbitals  u^,  Eq .  (20)  must  be  solved  iteratively, 
until  self-consistency  is  achieved.  Empirically,  if  one  takes  as  an  initial 
guess  that  p  is  approximately  a  sum  of  superimposed  atomic  charge  densities, 
then  the  convergence  of  this  procedure  is  fairly  rapid.  The  factor  of  2/3 
multiplying  the  potential  is  a  result  of  the  linear  dependence  of  Exc  on  p. 
This  also  has  a  consequence  that  the  Xa  eigenvalues  do  not  satisfy 
Koopman's  cheorem,  i.e.,  they  cannot  be  interpreted  as  ionization  energies. 
However,  it  can  be  shown  that  the  are  partial  derivatives  of  the  total 
expression  of  Eq.  (12)  with  respect  to  the  occupation  number, 


I 


de 

dni 


(21) 


If  E  were  a  linear  function  of  n^,  then  Koopman's  theorem  would  hold.  However, 
because  of  the  dominant  Coulomb  term,  E  is  better  approximated  by  a  quadratic 
function  in  n^.  This  leads  to  the  "transition  state"  approximation  which 
allows  one  to  equate  the  difference  in  total  energy  between  the  state  (n^,  n^) 
and  (n^-1,  n.+l)  to  the  difference  in  the  one-electron  energies  (Cj  -  e^)  cal¬ 
culated  in  the  state  (n^  -  1/2,  n^  +  1/2).  The  error  in  this  approximation  is 
proportional  to  third-order  derivatives  of  E  with  respect  to  n^  and  n j ,  which 
are  usually  small  (Ref.  31).  The  main  advantage  of  using  the  transition  state 
rather  than  directly  comparing  the  total  energy  values  is  computational  con¬ 
venience,  especially  if  -the  total  energies  are  large  numbers  and  the  difference 
is  small. 


The  relationship  of  Eq.  (21)  also  implies  the  existence  of  a  "Fermi  level" 
for  the  ground  state.  This  can  be  seen  by  varying  E  with  respect  to  n^  under 

the  condition  that  the  sum  In,  is  a  constant,  i.e., 

i  1 


(22) 
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This  implies  that  the 


implies  - —  =*  X,  where  X  is  a  Lagrangian  multiplier. 

total  energy  is  stationary  when  all  the  one-electron  energies  are  equal.  How¬ 
ever,  the  occupation  numbers  are  also  subject  to  the  restriction  0  *  n{  *  1. 
This  leads  to  the  following  conditions  on  the  ground  state  occupation  numbers 


(«i  <X0;  n.  =1 

(«i  >  X  •)*.  n,  =0 
«j  =  X  •)*  0<  Hj  <  | 


(23) 


In  other  words,  the  ground  state  eigenvalues  obey  Fermi  statistics  with 
X  representing  the  Fermi  energy.  It  should  be  noted  that,  in  contrast  to  the 
Hartree-Fock  theory,  where  all  the  n^  are  either  0  or  1 ,  the  Xa  model  predicts 
in  some  cases,  fractional  occupation  numbers  at  the  Fermi  level.  In  particu¬ 
lar,  this  will  occur  in  a  system  (such  as  transition  metal  or  actinide  atom) 
which  has  more  than  one  open  shell. 

The  Xa  model  differs  in  other  significant  ways  from  the  Hartree-Fock 
method.  In  fact,  the  simplification  introduced  in  approximating  the  total 
energy  expression  introduces  several  distinct  advantages  over  Hartree-Fock: 

1.  The  primary  advantage  is  purely  computational.  The  one-electron 
potential  in  Eq.  (20)  is  orbital- independent  anc  local,  i.e.,  it  is  the  same 
for  all  electrons  (except  in  the  spin-polarized  Xa  theory)  and  is  a  multipli¬ 
cative  operator.  On  the  other  hand,  the  Hartree-Fock  potential  is  nonlocal, 
or  equivalently,  there  is  a  different  local  potential  for  each  orbital.  This 
involves  a  great  deal  more  computational  effort,  especially  for  systems  des¬ 
cribed  by  a  large  number  of  orbitals.  It  has  been  shown  (Ref.  32)  that  the 
Xa  orbitals  for  the  first  and  second  row  atoms  are  at  least  as  accurate  as  a 
double-zeta  basis  set,  and  are  probably  better  for  larger  atoms  which  involve 
electrons  with  l  ^  2. 

2.  The  orbital-independent  Xa  potential  leads  to  a  better  one-electron 
description  of  electronic  excitations  of  a  system.  Both  the  unoccupied 

(n^  =  0)  and  occupied  (n^  =  1)  eigenfunctions  are  under  the  influence  of  the 
same  potential  resulting  from  the  other  N-l  electrons.  The  Hartree-Fock 
virtual  orbitals  see  a  potential  characteristic  of  the  N  occupied  orbitals, 
and  therefore  are  not  as  suitable  for  describing  the  excited  states.  Actually 
although  the  ground  state  virtual  eigenvalues  are  usually  a  good  description 
of  the  one-electron  excitations,  the  virtual  spectrum  of  the  transition  state 
potential  where  one-half  an  electron  has  been  removed  from  the  system  gives 
a  much  better  first-order  picture  of  these  levels  (Ref.  33). 


am 
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3.  As  has  been  shown  by  Slater  (Ref.  34),  the  Xa  model  rigorously 
satisfies  both  the  virial  and  He liman- Feynman  theorems,  independent  of  the 
value  of  the  parameter  a.  This  is  convenient  for  calculating  the  force  on  a 
nucleus  directly  in  terms  of  a  three-dimensional  integral,  rather  than  the 
six-dimensional  integrals  in  the  expression  for  the  total  energy  of  Eq.  (12). 
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7.  Computational  Aspects  of  the  Xa  Method 

In  application  of  the  Xa  model  to  finite  molecular  systems,  there  are  two 
practical  aspects  of  the  calculations  which  must  be  considered.  The  first 
concerns  the  choice  of  the  integration  framework  for  describing  the  molecular 
wavef unctions  and  the  second  deals  with  the  choice  of  the  exchange  parameter, 
a,  in  different  regions  of  space. 

In  computations  with  heteronuclear  molecules,  there  are  several  free 
parameters  that  must  be  chosen:  the  ratio  of  sphere  radii  for  the  atomic 
spheres  of  integration  at  a  given  internuclear  separation,  the  degree  of 
sphere  overlap,  and  the  value  of  the  exchange  parameter  in  the  atomic  spheres 
and  the  intersphere  region. 

It  has  been  found  that  changing  the  ratio  of  the  sphere  radii  for  the 
two  atoms  in  a  heteronuclear  diatomic  molecule  introduces  changes  in  the 
total  energy  that  can  be  large  on  a  chemical  scale  ('u  1  eV).  A  choice  for 
sphere  radii  based  on  covalent  bonding  radii  does  not  necessarily  provide  a 
good  estimate  for  these  calculations.  The  value  of  the  exchange  parameter, 
a,  and  the  sphere  radii  and/or  sphere  overlap  is  normally  fixed  in  Xa  cal¬ 
culations  for  crystals  where  the  geometry  is  fixed.  However,  to  develop 
a  potential  curve,  the  molecular  description  needs  to  change  substantially 
as  the  internuclear  separation  varies  and  the  changing  sphere  radii  include 
varying  fractions  of  the  total  molecular  charge  (Ref.  35).  Studies  made  at 
UTRC  have  shown  that  at  any  given  separation  the  total  energy  calculated  from 
the  Xa  model  is  a  minimum  at  the  radii  ratio  where  the  spherically  averaged 
potentials  from  the  two  atomic  centers  is  equal  at  the  sphere  radius. 


W  •  W 


(24 ) 


This  relationship  between  the  potential  match  at  the  sphere  boundary  and 
the  minimum  in  the  total  energy  appears  to  hold  exactly  for  "neutral"  atoms 
and  holds  well  for  ionic  molecular  constituents.  In  the  case  of  two  ionic 
species,  the  long  range  tail  of  the  potential  must  go  like  +2/R  from  one  ion 
and  -2/R  (in  Rydbergs)  for  the  other  ion  and  so  at  large  internuclear  separa¬ 
tions,  the  tails  of  the  potential  cannot  match  well.  However,  at  reasonable 
separations,  the  1/R  character  of  the  potential  does  not  invalidate  the  po¬ 
tential  match  criterion  for  radii  selection.  This  match  for  the  atomic  poten¬ 
tials  is  applied  to  the  self-consistent  potentials. 

In  molecules  with  significant  charge  sharing  in  the  bonds,  the  radii  of 
the  atomic  spheres  is  frequently  increased  in  Xa  calculations  so  that  an  over¬ 
lap  region  appears  in  the  vicinity  of  the  bond  (Ref.  36).  Studies  made  at 
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UTRC  show  that  the  contribution  to  the  total  molecular  energy  from  the  exchange 
integral  shows  a  minimum  at  the  optimum  sphere  radius  or  sphere  overlap.  This 
provides  a  sensitive  criterion  for  selecting  these  parameters. 


The  values  of  the  exchange  parameters  in  the  spherical  integration  region 
around  each  atomic  center  are  frequently  set  at  the  atomic  values  both  for 
neutral  and  for  ionic  molecular  constituents.  However,  for  light  atoms,  the 
value  of  a  which  best  reproduces  Hartree-Fock  results  varies  substantially 
with  ionicity.  In  argon,  the  following  table  compares,  for  the  neutral  atom 
and  the  positive  ion,  the  HF  energy  and  the  Xa  energy  calculated  for  several 
values  of  a. 


a 

Xu  Energy 

HF  Energy 

> 

O 

.72177 

526.8176 

526.8173 

Ar+1/2 

.72177 

526.5857 

_ 

.72213 

526.6007 

- 

Ar+1 

.72177 

526.2447 

.72213 

526.2596 

- 

.72249 

526.2745 

526.2743 

The  optimum  value  of  a  changes  even  more  rapidly  in  the  fluorine  atom,  going 
from  0.73732  for  F®  to  0.72991  for  F  Since  the  total  energy  depends 
linearly  on  u,  this  parameter  must  be  chosen  carefully. 


The  intersphere  exchange  coefficient  is  chosen  to  be  a  weighted  average 
of  the  atomic  exchange  parameters  from  the  two  constituents.  At  small  inter- 
nuclear  separations,  the  optimum  radius  for  an  atomic  sphere  frequently  places 
significant  amounts  of  charge  outside  that  atomic  sphere  -  charge  that  is 
still  strongly  associated  with  its  original  center  rather  than  being  trans¬ 
ferred  to  the  other  center  or  associated  with  the  molecular  binding  region. 

To  best  account  for  these  cases  the  weighting  coefficients  are  chosen  to 
reflect  the  origin  of  the  charge  in  the  intersphere  (or  outersphere  region). 


(Q„ 


<0 


O' 


-  o 


a. 


intersphere 


(«. 


0\ 


-  ) 


(25) 


where 

value 


U 

(Qs .  -  Qj  )  is  the  charge  lost  from  sphere  i  relative  to  its  atomic 
(or  ionic  value)  Q.°  and  a  is  the  atomic  exchange  parameter  for 
sphere  i.  This  value  for  ^intersphere  calculated  dynamically  -  it  is 
updated  after  each  iteration  in  the  self-consistent  calculation. 
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While  for  heavy  atoms,  these  changes  in  the  exchange  parameter  would  be 
small,  the  a  s  for  small  atoms  vary  rapidly  with  z  (and  with  ionicity).  The 
correct  choice  of  the  exchange  parameters  influences  not  only  the  total  energy 
calculated  for  the  molecule  but  also  in  some  cases  affects  the  distribution 
of  charge  between  the  atomic  spheres  and  the  intersphere  region. 


8.  Charge-Transfer  Calculations 

Even  though  low-energy  atom-atom  reactions  play  an  important  role  in  many 
physical  processes,  until  recently  comparatively  little  effort  has  been  devoted 
to  acquiring  a  knowledge  of  the  appropriate  cross-sections.  In  the  past, 
both  theoreticians  and  experimental  physicists  have  found  it  easier  to  study 
high-energy  collisions.  At  these  energies,  the  two  colliding  particles 
preserve  their  identities,  and  it  is  possible  to  treat  the  interaction  between 
them  as  a  perturbation.  There  is  no  guarantee  that  this  procedure,  known  as 
the  Born  approximation,  will  always  converge  to  the  correct  result  (Ref.  37). 

As  the  energy  of  the  colliding  particles  decreases,  it  is  necessary  to  take 
account  of  the  distortion  these  particles  undergo  during  the  collision.  The 
method  of  perturbed  stationary  states  was  developed  for  calculating  charge 
transfer  and  electronic  excitation  cross-sections  in  relatively  slow  collisions 
between  heavy  particles  (Ref.  38).  The  method  has  been  presented  in  both 
wave  and  impact  parameter  formalisms.  In  the  first  of  these,  the  entire 
system  is  treated  quantum-mechanically ,  while  in  the  latter  the  nuclei  are 
assumed  to  behave  as  classical  particles,  traveling  along  straight  line 
trajectories,  and  the  time-dependent  Schrodinger  equation  is  solved  to 
calculate  the  probability  of  various  types  of  electronic  transitions  (Refs.  39, 
40,  and  41).  Forcing  the  particles  to  travel  along  straight  lines  limits 
the  validity  of  the  impact  parameter  method  to  collisions  of  several  hundred 
electron-volts  or  greater  (Ref.  42).  The  wave  formulation  of  the  method  of 
perturbed  stationary  states  appears  to  be  one  practical  method  of  calculating 
thermal  energy  charge-exchange  cross-sections.  A  semi-classical  close-  ■*  ; \ng 
method  (Ref.  43),  based  on  an  average  scaling  procedure,  also  offers  ui'iity 
for  low  to  intermediate  collision  energies. 

The  use  of  the  wave  formalism  to  study  charge-transf er  reactions  at 
thermal  energies  dates  back  to  Massey  and  Smith's  pioneering  study  of  He+  +  He 
thermal  low-energy  scattering  (Ref.  44).  Strictly  speaking,  their  theory  is 
applicable  only  to  resonant  charge  transfer  reactions,  A  +  A+  -v  A+  +  A;  however, 
it  can  be  generalized  to  study  nonresonant  charge-transfer  reactions  as  well 
as  charge  transfer  into  excited  states.  If  the  origin  of  coordinates  is 
located  at  the  center  of  mass  of  the  nuclei  of  the  colliding  atoms,  the 
Schrodinger  equation  for  an  atomic  collision  can  be  written  in  center  of 
mass  coordinates  as 

(H-E)x  (F,R)  =  (h0-  v|  2  Vj-Vj-Ejx  (F,f?)-0  (26) 

where  r  represents  the  position  of  the  electrons,  R  is  the  vector  joining 
A  to  B,  Hq  is  the  Hamiltonian  for  the  system  when  the  nuclei  are  held  fixed, 

M  is  the  reduced  mass  of  the  two  nuclei.  Mis  the  sum  of  the  nuclear  masses, 
and  E  is  the  internal  energy  of  the  system,  including  the  electronic  energy. 
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Ignoring  heavy  particle  kinetic  energy  terms  in  the  center  of  mass  system 
results  in  a  modified  form  of  the  adiabatic  approximation  (Ref.  45)  and 
yields  perturbed  molecular  eigenfunction  \jjn(r,R)  which  satisfy  the  equation 


Here  en(R)  is  an  electronic  energy  level  of  the  molecule  perturbed  somewhat 
by  the  appearance  of  the  cross  teims.  The  wave  function  describing  the 
colliding  system,  X(r,R)  can  be  expanded  as 


x(?,R)*2  y/n(r,R)Fn(R) 


(28) 


The  various  scattering  cross  sections  are  determined  by  the  asymptotic  behavior 
of  the  Fn(R).  These  functions  are  determined  by  substituting  the  expansion 
Eq.  (28)  into  Eq.  (26).  Making  use  of  the  orthogonality  of  the  molecular 
eigenfunctions,  it  is  easy  to  derive  the  following  set  of  coupled  differential 
equations  for  the  Fn(R) 


(29) 


*2  <  ?»FmW} 

2  2 

where  V  (R)  and  k  /2M  are  the  potential  and  kinetic  energies  of  particles 
in  the  n  channel.  Many  of  the  difficulties  associated  with  trying  to 
calculate  thermal  energy  cross-sections  emanate  from  trying  to  derive  and 
solve  this  infinite  set  of  coupled  partial  differential  equations.  Until 
recently,  the  biggest  obstacle  in  the  calculation  of  low-energy  cross-sections 
has  been  the  inability  of  theorists  to  develop  accurate  molecular  eigen¬ 
functions.  For  those  problems  for  which  the  molecular  eigenfunctions  were 
available,  the  agreement  between  theory  and  experiment  has  been  very  reasonable 
for  the  amount  of  computational  effort  involved.  This  is  true  for  spin 
exchange  and  excitation  transfer  reactions  as  well  as  resonance  charge 
transfer  reactions  (Refs.  46-50). 
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Although  the  molecular  wave  functions  and  energies  available  in  the 
past  have  not  been  sufficiently  accurate  to  permit  extension  of  the  wave 
formalism  to  systems  having  more  than  about  four  electrons,  recent  advances 
in  calculational  techniques,  especially  for  two-center  systems,  have 
largely  overcome  this  problem.  In  particular,  recent  studies  have 
demonstrated  the  possibility  of  producing  highly  accurate  adiabatic 
electronic  wave  functions  for  systems  containing  as  many  as  40  electrons 
(Refs.  6,  17,  51,  52).  The  availability  of  these  small  but  flexible 
wave  functions,  which  have  the  property  of  connecting  with  the  correct 
separated  atomic  states,  increases  substantially  the  chances  for  successful 
and  practical  calculations  of  cross-sections  using  the  perturbed  stationary 
state  technique.  This  is  especially  true  of  charge  tranfer  into  excited 
states,  where  a  knowledge  of  a  number  of  the  low-lying  excited  states  is 
required. 

Even  with  a  reasonable  number  of  the  molecular  eigenfunctions,  the 
problem  of  calculating  cross  sections  is  far  from  solved;  this  is  especially 
true  of  charge  or  excitation  transfer  into  excited  states.  Many  of  the 
existing  studies  of  symmetric  resonance  reactions  are  based  on  the  two- 
state  approximation  which  limits  the  number  of  terms  in  the  expansion  Eq.  (14) 
to  two.  Under  these  conditions,  the  coupling  terms  on  the  right-hand  side 
of  Eq.  (29)  vanish  (Ref.  53).  We  have,  then,  only  to  solve  two  partial 
differential  equations,  instead  of  a  system  of  coupled  equations.  The 
situation  is  not  as  simple,  however,  for  nonsymmetric  reactions  or  for 
high-energy  collisions  involving  excitation.  For  these  problems,  the 
coupling  terms  are  the  source  of  the  transition  and  the  coupled  differential 
equations  have  to  be  solved  directly. 

While  a  great  deal  of  effort  has  gone  into  deriving  formal  theories 
of  inelastic  and  rearrangement  collisions,  relatively  little  work  has  gone 
into  trying  to  solve  the  resulting  equations.  This  is  especially  true  of 
low-energy  collisions,  where  the  lack  of  good  molecular  wave  functions  has 
prevented  one  from  evaluating  the  terms  coupling  the  different  channels. 
Previous  studies  which  were  part  of  the  UTRC  research  program  in  the 
electronic  structure  of  atoms  and  molecules  have  been  devoted  to  calculating 
matrix  elements  similar  to  some  of  the  terms  coupling  the  electronic  and 
nuclear  motion.  Some  of  the  required  work  involves  calculating  derivatives 
of  the  electronic  wave  functions  with  respect  to  the  internuclear  distance. 
This  task  is  made  simpler  by  the  use  of  compact  but  flexible  wave  functions 
such  as  those  studied  previously  at  this  center. 

The  biggest  obstacle  to  calculating  low-energy  cross  sections  is  the 
solution  of  the  infinite  system  or  coupled  partial  differential  equations 
describing  the  scattering,  Eq.  (29).  The  physics  of  the  problem  usually 
serves  as  a  guide  to  truncating  these  equations  to  a  system  of  finite  order. 

A  partial  wave  expansion  of  the  Fn(R)  leaves  a  large  number  of  sets  of 
coupled  differential  equations  (Refs.  54  and  55).  Since  the  number  of 
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equations  increases  as  the  collision  energy  increases,  there  is  no  single 
method  for  solving  these  equations.  At  thermal  energies,  a  direct  numerical 
integration  of  these  equations  is  feasible  (Refs.  56  and  57).  At  higher 
energies,  when  inelastic  collisions  and  charge  transfer  into  excited 
states  becomes  important,  the  trajectories  of  the  incident  and  scattered 
particles  may  be  nearly  classical  (Ref.  58).  Under  these  conditions,  it 
is  often  possible  to  use  the  WKB  wave  function  to  obtain  approximate 
solutions  to  these  equations  (Ref.  59). 
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9.  R-Matrix  Propagator  Solution  for  Charge  Transfer  Reactions 

The  theory  of  heavy  particle  charge  transfer  is  complicated  by  the 
fact  that,  in  general,  no  simple  natural  coordinate  system  exists  for 
which  a  uniform  description  of  reactant  and  product  states  can  be 
written.  Considering  the  reaction: 


0+[“*S]  +  N2  [X1!*]  -  NOMxV]  ♦  N  [ 4S) 


(30) 


The  internal  coordinates  describing  N,  cannot  be  uniformly  mass-scaled  to 
represent  NO  formed  as  product.  A  natural  coordinate  system  for  such  a 
collision  has  been  described  by  Marcus  (Ref.  60),  but  this  yields  very 
complicated  forms  for  the  kinetic  energy  operator  and  presents  enormous 
computational  difficulties. 

Recently,  Stechel,  et  al.  (Ref.  61)  have  described  an  approach  in 
which  two  separate  nonorthogonal  frameworks  are  employed  to  describe 
reactant  and  product  states.  The  total  scattering  wavefunction  is  written 
as: 


0  =  0a(R,Pj+  0y(R,Py) 


(31) 


where  a  and  y  represent  reactant  and  product  states,  respectively.  Expanding 
each  as  a  linear  combination  of  wave  functions  of  the  type: 


01  -  ^  f j  (R)  ) 

a 


(32) 


they  show  that  a  unique  connection  between  reactant  and  product  states  can 
be  achieved,  provided  the  basis  functions,  { 4>^  (p . )  }  are  taken  to  be  non¬ 
orthogonal,  i.e.,  Sa^  =  <<}>ot|<t>'<'>  i  0. 

We  note  that  Eq.  (32),  in  molecular  language,  is  a  product  wavefunction  of 
the  valence-bond  type.  We  also  note  that  an  SCF  or  MCSCF  framework  for  Eq.  (32) 
is  not  appropriate  for  treating  the  collisional  aspects  of  0+  +  N2  due  to  the 
orthogonality  constraints  that  must  be  imposed  on  such  wave funct ions .  The  use 
of  nonorthogonal  basis  functions  is  simplified  through  the  application  of 
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projection  algebra  techniques  for  the  construction  of  wavefunctions  of 
overall  symmetry.  Efficient  computer  codes  for  implementing  such  techniques 
has  been  developed  in  this  Research  Center  over  the  past  years  and  can  be 
utilized  in  constructing  the  proper  asymptotic  reactant  and  product  state 
channels. 

The  solution  of  the  charge  exchange  problem  can  be  accr.jplished  using 
a  direct  Numerov  integration  of  the  Schrodinger  equations  arising  from  the 
basis  representation  of  Eq.  (32).  Alternatively,  the  R-inatrix  propagator 
method,  originally  formulated  by  Diestler  (Ref.  62),  can  be  employed.  This 
approach  relates  the  known  asymptotic  solution  in  the  reactant  channel  to 
that  in  the  product  channel  through  a  series  of  recursive  steps.  Defining 
the  normal  derivative  to  the  surface  R  =  constant  by  Vnf(R,p),  we  seek  a 
propagator  operator  of  the  form: 


(33) 


The  explicit  matrix  representation  of  this  operator  can  be  written  as: 


tf^i  *  J.PJ 

*[Ri+^*P] 
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i  2  * 


r1  v  ib 
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(34) 


(33) 


Stechel,  et  al.  (Ref.  61),  have  given  explicit  expressions  for  constructing 
the  components  of  this  operator,  both  in  the  asymptotic  and  overlap  regions 
of  the  interaction.  The  initial  boundary  condition  in  the  reactant  channel 
is  R1  =  1/Xi,  where  Ai  is  the  eigenvalue  of  the  reactant  target  (in  this 
case,  N2) .  The  accuracy  in  developing  the  final  R-matrix  depends  on  the  chosen 
step  size,  h_^,  which  is  taken  as  an  adjustable  parameter  that  can  be  optimized 
during  the  course  of  the  propagator  solution. 

We  have  recently  implemented  this  R-propagator  analysis  at  UTRC,  using 
a  modified  version  of  the  code  developed  by  Stechel,  et  al.  (Ref.  61).  One 
application  has  been  the  low-energy  collisional  de-excitation  reaction: 


N*[  2D]  +  0(3P]  -  NI4S)  +  0[3P] 


(36) 
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which  occurs  along  a  curve-crossing  21+  -  2£+  potential  surface.  The  results 
of  this  study  have  been  reported  in  AFGL-TR-79-0190  (Ref.  12)  and  are  presently 
being  prepared  for  journal  publication. 


10. 


Spin-Projected  Unrestricted 


Hartreo-Fock  Method 


The  unrestricted  Hartree-Fock  (UHF)  method  developed  by  Pople  and  Nesbet 
(Refs.  63  and  64)  yields  the  best  single-determinant  approximation  to  the  exact 
wavefunction  for  an  atomic  or  molecular  system.  Such  a  wavefunction  incorporates 
correlation  by  allowing  orbitals  of  different  spin  to  adjust  to  spatially  different 
forms,  thus  breaking  the  symmetry  restrictions  of  the  conventional  (RHP)  method 
(Refs.  15  and  16).  It  is  necessary,  however,  to  project  from  sucli  a  wavefunction, 
a  properly  antisymmetrized  spin  and  angular  momentum  state  in  order  to  define 
eigenstates  and  e igenenergies  corresponding  to  observable  spectroscopic  states. 


Let  ja>,  1 6>  refer  to  the  doubly-occupied  molecular  orbitals  (MO's),  'P  » 
4>g.  Let  Iy>,  1 6 > ,  j  t >  refer  to  the  singly  occupied  MO's.  We  assume  that  all 
MO's  have  been  subjected  to  a  transformation  to  orthogonal  form.  Let  CL.^  be 
the  coefficient  associated  with  the  permutation  y  <->■  6  (this  permutation  is 
(-1)  times  the  overlap  of  the  permuted  spin  eigenfunction  with  the  original 
spin  eigenfunction),  and  adopt  the  convection  =  -1.  Then  the  expectation 
value  of  the  nonrelativistic  Hamiltonian 

3?--  l  U(i)  +  Z  V(.,j)  (37) 

1  '<) 


is  given  by 


<3f>  =  2  Z  <a|u|a  >  +  Z  <  y  |u|  r  >  +  2  Z  <  a/3  |v|  a/?> 

<a^|v)/Sa>  +  2  Z  <ar|v|ar>  - 1  <a  y  Ivl  ya> 

<*P  ay  1  ay  1  1  (38) 

fT^8<rS  lvl  \  Zg  Cn  <X8  |v|  Sy> 


where  i,  j  refer  to  electron  numbers,  a,  6..  y,  6  etc.,  refer  to  MO's  and  the 
sums  run  over  all  a,  B.... 
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The  equations  for  determining  the  optimum  MO's  can  be  derived  from  the 
following  variational  form 


8 


«ya<a 


\Y>  +  €, 


ay 


</|a>)- 


I 

ys 


ys 


y|8: 


=  0  (39) 


These  equations  are 


(U+2  Jc-Kc  +  J„-  ±K0)|a>  =  |  +  i|%|8>  <10> 

(u  +  2  jc-Kc  +  Jo)|y>  +  |cr8|<s|v|y>s>  =  |‘flyl^>  +  |vls>  <41) 

where 

|<8|v|y>8>  =  /dV$*  (7)  vCT,?)^>y(T)<^(7) 

<w|v|w'>  =  f  d*'  w+(T)  v(T.T')  w'(T) 

Jc  =  |  </3|v|/9>  J0=|<8|v|8>  (42) 

kc|w>=2  |  </3|v|w>  (3  >  k0|w>=|  |<  8 1  v|w>  8  > 

Left-multiplying  Eq.  (4o)  by  <y\  and  Eq .  (41)  by  <u|,  these  equations  can  be 
combined  to  yield 

*ya  =  |{2Cyg  +  i)  <ys|v|8a>  (43) 

Noting  that 


28 


Finally,  Eqs.  (48)  and  (49)  can  be  rearranged  to  yield  the  same  Hermitian 
operator  on  the  lbs  with  the  resulting  definition  of  the  one-electron  eigen- 
energies  for  the  closed  and  open-shell  eigenstates: 


«a  =  «*|u  +  2Jc-  Kc  +  J0-  ^K0la> 


(50) 


fy=  <y|u  +  2JC-  Kc  +  JQ-  ^K0+  -^M+|y> 


(51) 


Equations  (50)  and  (51)  represent  the  correct  one-electron  energy  expressions 
for  the  spin-projected  eigenstate. 

This  formalism  has  been  developed  to  extend  the  conventional  HF  method 
to  include  split-shell  correlation  and  proper  spin  and  symmetry  projections. 
It  has  been  incorporated  into  a  computer  program  using  Gaussian-type  orbitals 
(GTO's)  as  the  elementary  basis  functions.  A  crucial  feature  of  this  method 
is  that  dissociation  always  follows  the  lowest  energy  pathway  thereby  permit¬ 
ting  a  proper  description  of  bond  formation  and  breakage.  In  contrast,  RHF 
or  MC-SCF  methods  often  exhibit  improper  dissociation  character  or  exhibit 
size  inconsistencies  owing  to  correlation  energy  changes  in  going  from  molec¬ 
ular  geometries  to  separated  atom-molecule  or  atom-atom  dimensions. 


DISCUSSION  OF  RKSUI/l’S 


Background 


For  several  years,  this  Center  has  been  involved  in  detailed  ab  initio 
quantum  mechanical  studies  of  the  electronic  structure,  radiative  and 
kinetic  properties  of  diatomic  and  polyatomic  molecules.  These  studies 
have  focused  on  air  molecules  and  metal  oxides,  with  applications  to  weapons 
technology. 

Our  initial  quantum  mechanical  studies  dealt  with  the  electronic 
structure  and  radiative  lifetimes  of  the  air  system  and  certain  diatomic 
metal  oxide  species  which  may  be  of  abundance  in  the  perturbed  atmosphere 
following  a  thermonuclear  explosion.  These  studies  have  included  the  systems 
NO,  N0+,  NO2 ,  N2O,  TiO,  and  TiO+.  Such  species  commonly  exhibit  strong  LW1R 
radiation  and  a  quantitative  measure  of  this  effect  is  needed  to  analyze 
the  atmospheric  kinetics.  These  systems  are  of  direct  relevance  to  current 
DNA  interests  and  to  the  general  problem  of  identifying  atmospheric  kinetics 
and  radiation  properties.  The  studies  were  carried  out  using  methods  of 
molecular  quantum  mechanics  coupled  with  available  spectroscopic  information 
on  excited  electronic  states. 

As  a  result  of  these  investigations,  which  were  mainly  carried  out  under 
Air  Force  sponsorship,  this  Center  has  available  a  large  data  base  of 
calculated  electronic  energies  and  transition  probabilities  for  air  system 
molecules  (and  molecular  ions)  and  for  metal  oxide  systems.  This  base 
includes  calculations  for  the  ground  and  many  electronically  excited  states 
of  these  molecules.  Several  of  these  systems  have  been  analyzed  in  detail 
(Refs.  65-75).  To  the  best  of  our  knowledge,  no  other  laboratory  in  the 
country  has  such  an  extensive  catalogue  of  theoretical  calculations  of 
the  electronic  structure  of  the  air  system  and  of  the  light  metal  oxides  or 
our  computational  experience  with  these  types  of  molecules.  In  addition, 

UTRC  has  fully  developed  sophisticated  computer  programs  for  studying  both 
diatomic  and  polyatomic  molecules.  These  programs  have  been  extensively 
used  at  UTRC  and  at  the  Air  Force  computational  laboratories  at  AFWL, 

Kirtland  Air  Force  Base,  and  AFGL,  Hanscom  Field. 

As  a  continuation  of  our  theoretical  research  on  atmospheric  chemistry 
under  the  subject  contract  F19628-80-C-0209 ,  we  have  examined  the  reaction 
kinetics  for  the  ion-molecule  reaction  0+  +  ^  -*■  N0+  +  N.  /  articular 
reaction  which  is  very  important  in  N0+  formation  is: 

o+[*s]  +  Na*[x'z;]-TO*[x'3:+] ,  n[*s]  (52, 

AE  =  -  1 1  ev 
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Another  ion-molecule  reaction  whose  role  is  less  certain  is: 


°14s]  +  n*[a3i;]~no+[x'Z+]  +  n[4s] 

AE  =  -73  ev 


(53) 


Analysis  of  the  kinetics  of  these  reactions  requires  the  calculation  of  the 
appropriate  reaction  surface  or  surfaces  followed  by  a  kinematic  analysis 
of  the  motion  of  the  reactants  on  such  surfaces.  In  this  work,  we  are 
primarily  concerned  with  analysis  of  the  0+  +  Nlf  reaction  (l.q.  52)  in  an 
energy  region  0.1  -+  10.0  eV,  with  this  energy  either  as  center-of-mass  (CM) 
translational  energy  or  as  vibrational  energy  of  the  N2  molecule.  Our  first 
concern,  therefore,  is  to  seek  possible  reaction  paths  which  strongly 
couple  reactants  to  products.  We  first  examined  states  of  quartet  symmetry 
which  exhibit  allowed  connections. 

Potential  Energy  Surfaces 

The  molecular  correlations  that  connect  reactant  to  product  states  are 
given  in  Tables  1  and  2.  The  correlation  diagrams  shown  in  Tables  1  and 

2  list  the  reactant  and  product  states  in  increasing  energy  and  show  the 
symmetry  of  the  electronic  states  of  for  unilinear  collision  only.  The 

lowest  quartet  state  in  a  linear  (C  )  0+  -  ^  conformation  corresponds  to 
symmetry.  In  Figs.  1  -  3,  we  show  the  adiabatic  correlations  for  the  quartet 
states  of  ^0+  in  all  possible  symmetries,  linear  (Crov)  ,  bent  (Cg) ,  and 
symmetric  perpendicular  (C2V) •  In  a  bent  geometry,  the  lowest  potential 
energy  surface  has  ^A"  symmetry  but  details  of  the  energetics  of  this  surface, 
to  include  possible  avoided  crossings  with  higher  states,  require  accurate 
quantum  mechanical  calculations  of  the  N20+  molecular  ion  as  a  function  of 
geometry.  A  limited  region  of  this  lowest  potential  energy  surface  can  be 
described  using  a  multipole  expansion  analysis  for  the  long-range  interactions 
as  described  below. 


Long-Range  1 nterac tions 


The  long-range  energy  of  interaction  for  this  system  can  be  conveniently 
represented  in  the  natural  coordinates  for  the  collision  as  illustrated 
in  Fig.  4.  The  collision  angle  6  can  be  related  to  the  molecular  angle  y 
through  use  of  the  relations: 


[ron2  +  (-^rf -R0nrnn  C0Sy] 

sm  (3 


% 


sm 

x 


(54) 
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In  terms  of  x  and  3,  the  long-range  energy  of  interaction  of  0+  +  N2 
can  be  written  as  (Refs.  76  and  77): 


V  (0*  +  N?  )  ~  +  Q2  z  COS2/3-l) 


(a||(N3)COS2  /3-t-q  j_(N2)  Sin2#) 


3a(04)(q||{N2)cos2/3  +  a  j_(n2)  sin2/?) 
2  x6  (I‘,(0<)+  I ' '( N2) ) 


The  first  term  in  Eq.  (55)  represents  the  ion-quadrupole  interaction,  the 
second  term  represents  the  main  component  of  the  energy  of  induction 
between  the  charge  and  the  induced  dipole  and  the  third  term  is  the  London 
dispersion  energy  based  on  the  approximate  Drude  formula.  A  small  con¬ 
tribution  (x-8)  to  the  energy  of  induction  has  been  neglected  owing  to 
the  smallness  of  the  permanent  moments  for  this  system. 

Tn  the  product  channel,  the  long-range  energy  of  interaction  of 
N  +  N0+  can  be  written  as 


V  ( N  +  NO* )  « 

X—  00 


3a(N)(a  ||(N04)  cos2/?  +  ax(NOMsin2/?) 
2  x6(  I'1  (N)  +  r'(NCn) 


where  3  refers  to  the  corresponding  collisional  angle  measured  in  the 
product  channel  framework.  Note  that  the  ion-quadrupole  interaction 
vanishes  here  since  the  N  atom  is  in  a  state. 

In  Table  3,  we  give  values  for  the  parameters  appearing  in  Eqs.  (55)  and 
(56).  Since  dipole  polarizabilities  for  the  ions  are  not  available,  they 
are  estimated  from  the  corresponding  neutral  species  using  expectation 
values,  r?,  calculated  from  corresponding  neutral  and  ionic  Hartree-Fock 
wave functions . 


Values  of  the  long-range  interactions  calculated  from  Eqs.  (55)  and 
(56)  are  presented  in  Tables  4  and  5,  respectively.  Comparison  with 
interaction  potentials  derived  from  an  MC-SCF  calculation  (Ref.  78) 
indicates  significant  differences.  This  suggests  an  inaccurate  representation 
of  polarization  forces  using  this  MC-SCF  wavef unction . 

Ab  Initio  Calculations 


The  first  ab  initio  calculations  for  the  ^0"’’  system  are  those 
reported  by  Pipano  and  Kaufman  (Refs.  79  and  80)  for  a  col  linear  arrangement 
of  the  nuclei.  These  results  were  obtained  using  an  atom-optimized  minimum 
Slater-type  orbital  (STO)  basis  with  a  very  large  configuration-interaction 
(Cl)  constructed  by  including  all  single  and  double  electron  excitations 
from  the  ground  reference  state.  The  final  Cl  sizes  (1379  for  ^j]  symmetry  and 
1310  for  t*'i~  symmetry)  represented  pioneering  work  for  a  molecular  system  of 
this  complexity.  Pipano  and  Kaufman  found  that  both  the  **T>  and  channels 

lie  too  high  in  energy  to  permit  a  direct  reaction  mechanism  from  0+  +  N2  to 
NO+  +  N.  Instead,  they  postulated  that  this  reaction  must  involve  multi- 
surface  mechanisms  involving  the  crossings :  <■  -»  or  possibly  ^ l~  -*• 

As  we  shall  see  below,  for  low  collisional  energies  (<  0.1  eV ) , 
the  multi-surface  reaction  of  Pipano  and  Kaufmann  is  the  most  probable 
mechanism  but  the  experimental  evidence  supports  a  direct  reaction  mechanism 
for  higher  collisional  energies. 

A  limited  study  of  bent  conformations  of  N'20+  was  also  reported  by 
Pipano  and  Kaufman  but  an  unfortunate  choice  of  the  bending  angle  (5° 
relative  to  a  linear  conformation)  led  them,  as  we  shall  see,  to  the 
incorrect  conclusion  that  the  energy  increased  with  bending  the  molecule. 

F.xtensive  ab  initio  SCF  and  MC-SCF  studies  of  the  lowest  potential 
surface  of  N90+  have  been  reported  by  Hopper  (Refs.  81  ar.d  82).  Hopper 
found  that  a  bent  conformation  of  the  nuclei  for  ^O-1"  resulted  in  a 
significant  lowering  of  the  energy  from  the  collinear  arrangement  studied 
by  Pipano  and  Kaufman.  He  reports  a  saddle  point,  with  a  barrier  height  of 
2.6  eV,  in  the  lowest  surface  at  the  geometry:  RN^  =  2.A2  bohrs,  R^q  = 

2.38  bohrs,  c.  NNO  =  120°.  These  bond  lengths  represent  an  expansion  from 
the  equilibrium  values  for  N0+  (2.0092  bohrs)  and  N2  (2.0742  bohrs).  In 
another  study  (Ref.  78),  Hopper  reports  calculations  on  the  lowest  ^YT 
surface  resulting  from  a  collinear  ariangment  of  the  nuclei.  Again,  he 
finds  a  saddle  region  at  considerably  expanded  bond  lengths  (R^  =  2.8  bohrs, 

=  2.6  bohrs)  and  a  calculated  barrier  height  of  't  8.0  eV  relative  to 
the  0+  +  k2  asymptote. 

Both  of  these  surfaces  exhibit  barrier  heights  that  are  much  too 
large  to  explain  the  experimentally  known  energy  dependence  of  the 
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0+  N,  -*•  N0+  +  N  reaction.  Hopper  argues  (Ref.  82)  that  a  complicated 
reaction  coordinate  begining  with  a  polarization  ^Z  state  and  switching 
to  a  bent  conformation  at  the  saddle  point  yields  a  reaction  path  that 
has  a  barrier  less  than  'v  0.2  eV.  Such  a  value  would  be  in  accord  with 
the  best  experimental  evidence.  His  argument  for  this  low  barrier  height, 
however,  is  not  based  on  his  calculated  value  (2,6  eV)  but  rather  on  a 
circumlocutory  thermodynamic  cycle  argument  involving  a  mix  ■  I  calculated, 
experimental  and  estimated  data.  It  was  felt,  therefore,  that  a  re¬ 
examination  of  the  lowest  potential  surface  for  N20+  should  be  carried  out. 

Open- She ^n^tiO-Cal^u^ations 

Our  initial  studies  were  directed  toward  an  examination  of  the  lowest 
^1“  and  H  surfaces  of  N20+  using  a  large  Cl  and  employing  accurate  Slater- 
type  orbit  (STO)  basis  sets.  This  Center  has  operational  a  very  sophisticated 
computer  code  for  polyatomic  calculations  which  incorporates  these  features. 
Boi.4  single-zeta  (minimum  basis)  and  contracted  double-zeta  STO  basis  were 
examined.  In  either  case,  a  complete  Cl  represents  1488  configurations  for 
and  531  configurations  for  These  Cl  expansions  include  all  possible 

orbital  occupancies  and  all  possible  spin  coupling  among  the  valence 
electrons . 

Previous  studies  in  this  Center  on  N20+  (Ref.  83)  had  shown  that  SCF 
or  MC-SCF  calculations  may  exhibit  artifically  high  reaction  barriers  in 
C<„v  symmetry.  The  origin  of  the  problem  is  that  SCF  or  MC-SCF  methods 
generally  employ  maximum  pairing  in  their  MO  structure  and  thus  impose 
constraints  on  spin-recoupling  that  may  occur  along  the  reaction  coordinate. 
The  SCF  structure  of  N2o+[^I-]  can  be  written  as 


,la^  2a^  30^40^50^6  o^l  ^7  8  < 


where  at  short  0  -  N2  separations,  corresponding  to  product  ion  formation,  we 

can  identify  the  principal  AO  composition  of  these  MO's  as  follows: 


la,  2a,  3o  = 
4o  = 
5o  = 
6o  = 
7a  = 
8o  = 
9o  = 

ITT  = 
2  tt  = 
371  * 


2s0 

2sfj(outer) 
2sjg(inner) 
2po (N0+) 
2poN(outer) 
2p*  (N0+) 

2px  (N0+) 
2pTTw  (inner) 
2p*  (N0+) 


At  large  O’*"  -  N2  separations,  the  MO  composition  for  the  2p  electrons  must  be: 

7o  ■=  2po(N2) 

80  «=  2po(0+)  (59) 

lx  =  2px(No) 

2x  =  2pTt(0+) 

Thus  both  the  a  and  it  orbitals  are  completely  recoupled  spin-wise  in  going 
from  0+  +  N£  to  N  +  NO  .  In  Fig.  5,  we  present  a  general  vector  coupling 
diagram  which  indicates  that  there  are  48  separate  spin-states  for  nine 
electron  quartet  couplings.  The  reaction  0+  +  N2  -*■  N0+  +  N  involves  a  complete 
spin-recoupling  as  the  pairing  changes  in  going  from  reactants  to  products, 
which  suggests  that  a  valence  bond  treatment  is  more  appropriate  than  SCF 
or  limited  MC-SCF  expansions.  Previous  studies  on  the  much  simpler  surface 
have  indicated  the  importance  of  including  higher  spin-coupled  states 
(Ref.  84). 

Preliminary  studies  were  carried  out  under  AFGL  Contract  F19628-77-C- 
0248  which  confirm  that  the  spin-recoupling  that  occurs  in  going  from 
reactants  to  products  is  not  adequately  represented  within  SCF  or  limited 
MC-SCF  frameworks.  A  calculation  for  the  lowest  state  of  N20+  at 
intermediate  internuclear  separations  yields  a  spin-correlation  error  of 
6  eV.  This  value  closely  corresponds  to  the  8  eV  barrier  found  in  the 
MC-SCF  calculations  by  Hopper  (Ref.  78).  The  major  error  arises  from  the 
7 o  and  lir  molecular  orbitals  which  change  in  composition  from  primarily 
N2  in  character  at  large  0+  +  N2  separations  (reactants) ,  to  mainly  N0+ 
at  large  N  +  NO+  separations.  This  error  corresponds  to  2  eV/ electron 
pair,  a  value  that  is  typical  of  many  atomic  and  molecular  systems. 

A  second  problem  with  limited  extent  MC-SCF  calculations  is  size 
inconsistency  errors  arising  from  changes  in  the  correlation  energy  in 
going  from  N20+  to  0+  +  Nj  (N0+  +  N)  and  to  the  (0+  +  N  +  N)  separated 
atom  limits.  The  SCF  wavef unction  given  in  Eq.  (57)  dissociates  properly 
in  both  the  0+  +  N2  and  N0+  +  N  limits  but  is  incapable  of  proper 
description  of  the  plateau  dissociation  to  0+  +  N  +  N.  In  this  region, 
the  minimum  Cl  (or  MC-SCF)  must  include: 
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The  first  configuration  is  the  basic  SCF  representation  and  the  others  are 
those  required  for  correct  atomic  limit  behavior.  Hopper  (Ref.  78)  has 
included  configurations  (1),  (2),  (3),  and  (5)  [his  (1),  (2),  (3),  and 
(13)]  in  his  MC-SCF  construction  for  the  lowest  state  but  these  are 
insufficient  for  proper  atomic  dissociation.  As  a  result,  it  is 
predictable  that  his  MC-SCF  calculations  will  exhibit  increasing  correlation 
errors  for  N20+  geometries  tending  toward  Qfr  +  N  +  N  dissociation.  Thus  the 
saddle  region  (extended  N-N  and  N-0  bonds)  will  be  poorly  described  relative 
to  0+  +  N2  (N0+  +  N)  using  Hopper's  MC-SCF  construction.  This  is  borne  out 
by  his  calculated  energy  for  the  saddle  point  on  the  surface.  His 
reported  energy  at  the  barrier  is  -183.160  hartrees,  a  value  already  equal 
to  the  atomic  dissociation  limit  (-183.159998)  with  his  basis,  set.  Hopper's 
predicted  location  and  energetics  of  the  barrier  region  must  therefore 
reflect  this  correlation  dissociation  error. 

The  results  of  our  CI-STO  basis  calculations  for  symmetry  predict 
a  saddle  point  at  Rj^  ^  R^0  a,  2.2  bohrs  with  a  barrier  height  of  0.1363 
hartrees  (3.7  eV)  relative  to  the  asymptote  of  0+  +  N2>  This  Cl  consisted 
of  a  full  expansion  of  all  a  and  u  valence  space  of  N20+  and  included  all 
eight  configurations  listed  above.  Based  on  this  result,  we  concluded  that 
a  linear  conformation  of  nuclei  is  not  the  lowest  possible  quartet  surface 
for  N20+. 

A  similar  CI-STO  basis  calculation  was  carred  out  for  symmetry. 

At  «  RqN  *  2.0  bohrs,  the  lowest  ^fl  state  lies  0.1509  hartrees  above  the 
0+  +  N2  asymptote,  very  close  to  the  ^ L~  surface.  This  state  has  a  dominant 
configuration  that  can  be  represented  as 
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[la2  ...  6 o2  7c2  In3  2tt3  3tt] 


(60) 


At  shorter  0  -  N2  separations,  the  dominant  configuration  switches  to 

4n  [lo2  ...  6o2  7 o2  lit4  2ti 2  3tt ]  (61) 

and  at  large  0+  -  N2  separations,  the  dominant  configuration  is 


4II  [lo2  ...  6a2  7o  8a  In4  2n3]  .  (62) 

Thus  a  complicated  crossing  of  the  surfaces  occurs  between  R^_q  =  2.0  - 
3.0  bohrs.  However,  none  of  these  linear  surfaces  is  low  enough  to  explain 
the  observed  energy  dependence  of  this  reaction. 

The  results  of  these  linear  conformation  studies  led  us  to  examine 
the  lowest  4A"  surface  at  several  internuclear  separations  and  bond  angles. 
Owing  to  cost  and  time  considerations,  the  ^A"  surface  was  studied  in  a 
projected  UHF  framework  using  Gaussian-type  orbitals  (GTO's)  in  single-?, 
double-?  and  double-?  +  single  d-polarization  function  approximations.  The 
latter  two  expansions  were  checks  on  the  work  of  Hopper  (Refs.  78  and  82). 
After  much  experimentation,  it  became  clear  that  the  minimum  energy  reaction 
path  lies  close  to  a  45°  trajectory  of  0+  +  N2  as  illustrated  in  Fig.  4.  Of 
course,  all  possible  trajectories  from  X  =  0  -+  X=90°  will  be  traversed,  but 
only  those  near  the  45°  channel  lead  to  the  reactive  rearrangement  collision: 
0+  +  N2  ->  N0+  +  N . 

To  insure  a  systematic  comparison,  the  linear  surfaces  for  the  lowest 
l*y~  and  4j{  states  were  recalculated  within  the  UHF  approximation.  They  are 
given  in  Table  6  and  illustrated  in  Fig.  6.  Our  calculated  barrier  for  the 
lowest  4"“  surface  is  0.128  hartrees  (3.5  eV) ,  close  to  the  value  found  for 
the  CI-STO  basis  calculations.  A  lower  linear  reaction  pathway  involves  a 
]1  -  4i“  surface  crossing  at  an  energy  of  'v  0.11  hartrees  (3.0  eV) ,  a  value 
still  too  high  to  explain  the  collisional  data.  A  comparison  of  the  several 
ab  initio  calculations  of  the  lowest  4£~  surface  is  shown  in  Fig.  7.  We 
have  also  included  the  analytical  long-range  potential  for  comparison  with 
the  previous  studies  of  Pipano  and  Kaufman  (Ref.  80)  and  Hopper  (Ref.  78).  In 
general,  Pipano  and  Kaufman's  CI-STO  results  tend  to  over-emphasize  the  long- 
range  behavior  whereas  Hopper's  MC-SCF  results  exhibit  an  unbalance  in  their 
long-range  behavior.  The  excess  long-range  attraction  exhibited  by  Hopper V 
calculation  in  the  0  +  N9  channel  is  most  probably  a  reflection  of  an 

unbalance  in  errors  in  his  calculations  for  N2  [X-*-Ig  ]  (2.03  eV  error  in  Dg) 
and  NO+  [x^T+]  (3.79  eV  error  in  De) .  The  extra  attraction  in  the  0+  +  N2 
reactant  channel  and  the  extra  repulsion  in  the  N0+  +  N  product  channel  are 
a  reflection  of  these  asymptotic  errors  in  going  from  the  reactants  0+  +  N2 
to  N20  and  finally  to  products  N0+  +  N.  However,  all  calculations  clearly 
illustrate  that  linear  surfaces  cannot  explain  the  energy  dependence 
of  this  reaction. 
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As  mentioned  above,  our  lowest  energy  reaction  path  closely  follows  a 
45°  trajectory.  An  optimization  of  the  saddle  region  along  this  path  was 
carried  out.  This  reaction  surface  is  in  a  convenient  coordinate  frame¬ 
work  for  collisional  studies  but  differs  from  conventional  studies  which 
use  the  bond  angle  of  the  molecule  as  an  independent  surface  geometry 
parameter.  Our  constant  trajectory  angle  of  45°  corresponds  to  a  cut 
through  many  surfaces  with  the  molecular  bond  angle  changing  from  about 
115  -*  120°.  The  lowest  energy  reaction  path  is  given  in  Table  7  and  is 
illustrated  in  Fig.  8.  We  find  a  barrier  height  of  n,  0.01  hartrees  (0.3  eV) 
and  a  relatively  flat  surface  in  the  region:  3.4  <  x  <  3.8  bohrs.  Also 
shown  in  Fig.  8  is  the  reaction  path  calculated  by  Hopper  (Ref.  82)  within 
a  restricted  HF-SCF  framework  and  at  a  constant  molecular  angle  of  120° . 

This  translates  to  a  collisional  angle  varying  from  about  40  -►  45° ,  close 
to  our  optimum  trajectory.  Hopper's  basic  conclusion  that  the  lowest 
energy  quartet  surface  for  N20+  is  bent  is  thus  confirmed  by  these  studies. 
The  0+  +  Nt2  -*  N0+  +  N  reaction  would  appear  to  take  place  on  a  single 
hypersurface  for  collisional  energies  between  “v  0.1  and  several  electron 
volts,  namely  for  collisional  energies  sufficient  to  cross  the  barrier  in  the 
lowest  "A"  surface. 

Potential  Surface  Modeling 

There  are  several  approaches  that  can  be  taken  to  yield  a  functional 
representation  of  the  reaction  potential  energy  surface.  The  ab  initio 
calculations  produce  a  large  table  of  numbers  that  must  somehow  be  fitted 
and  smoothed  for  use  in  trajectory  calculations  of  reaction  rates.  Our 
first  approach  was  to  examine  an  expansion  of  the  interaction  potential 
between  0+  and  N2  in  terms  of  Legendre  polynominals  as 


V(R,x,  X)  -  2  Vj  (R,  x)  F|  (cos  X)  (63) 

even 

This  form  has  been  examined  previously  for  Li+  +  H2  (Ref.  85)  and  H“  +  H2 
(Ref.  86)  collisional  studies.  An  analysis  of  the  lowest  ^A"  surface  of 
N20+  in  terms  of  Eq.  (63)  yields  a  slowly  convergent  expansion  owing  to 
complications  of  representing  both  the  attractive  long-range  potential 
for  linear  0+  +  N2  and  the  relatively  flat  saddle  region  for  >  =  45°.  Thus 
the  utility  of  employing  Eq.  (63)  rather  than  a  direct  table  interpolation 
scheme  seems  questionable. 

An  alternate  approach  is  to  choose  a  more  restricted  functional  form 
than  that  represented  by  Eq .  (63)  but  one  that  can  easily  incorporate  the 
known  asymptotic  behavior  of  the  system  at  the  0+  +  N2  and  rtO+  +  N  limits. 

The  London-Eyring- Polany i-Sato  (LEPS)  model  (Ref.  87)  has  been  very  successful 
in  applications  to  atom-molecule  collisions  (Ref.  88).  For  the  N20+  system, 
the  three-body  interaction  potential  can  be  written  as 


vabc  '  P  [^ab  +  ^bc  +  ^ca] 


(64) 

~  yf  [^ob  "  abc^  +  (Qbc~  aoc^2  +  (Qac  "  aob^2] 


where  and  are  Che  Coulomb  and  exchange  two-body  interactions  and  p 
is  treated  as  an  empirical  parameter.  The  Q  —  and  a. ,  are  related  to  the 
diatomic  potential  functions  by  1'* 


V,j  =  °ij  +  ^a,, 


(65) 


where  the  V. .  can  be  expressed  by  some  empirical  form  such  as  a  Morse 
function  (Re^ .  89).  For  these  present  studies,  Eq.  (65)  was  fitted  to 
a  Hulbert-Hirschfelder  potential  function  (Ref.  90)  modified  to  include 
the  long-range  polarization  tail.  The  adjustable  parameters  were  used 
to  fit  Eq.  (64)  in  critieaL  regions  of  the  potential  surface  such  as 
the  barrier  and  asymptotic  regions.  The  fitted  minimum  energy  LEPS 
reaction  path  is  close  to  the  UHF  surface  in  the  region  of  the  saddle  point 
as  can  be  seen  in  Fig.  8. 

In  Figs.  9  and  10,  we  give  the  LEPS  representation  for  the  ^O"1*  surface 
in  C2v  symmetry.  Figure  10  illustrates  the  large  barrier  for  reaction  in 
this  geometry.  The  collinear  LEPS  surface  is  shown  in  Fig.  11  which  again 
illustrates  a  large  reaction  barrier.  A  LEPS  surface  cut  at  a  fixed  molecular 
angle  of  120°  is  shown  in  Figs.  12  and  13.  The  barrier  height  is  much 
reduced  for  this  bent  geometry.  Finally,  in  Fig.  14,  we  show  the  minimum 
energy  reaction  surface  that  was  employed  in  the  collisional  studies.  An 
important  feature  of  this  surface  is  that  the  saddle  point  is  nearly 
symmetrically  located  between  the  entrance  and  exit  channels.  Thus  either 
translational  or  vibrational  energy  should  be  effective  in  this  reaction. 
Translation  energy  is  effective  in  surmounting  the  barrier  while  vibrational 
energy  is  effective  is  converting  to  N0+  +  N  products  in  the  exit  channel. 

A  barrier  region  deep  in  either  the  entrance  or  exit  channels  would 
exhibit  bias  toward  translational  or  vibrational  energy,  respectively.  The 
surface  represented  by  Fig.  14  is  that  employed  in  the  collisional  studies 
described  below. 

Kinetics  of  the  0+  +  N->  Reaction 

Our  kinetic  studies  were  restricted  to  the  analysis  of  the  reaction 

0*[*s]  +  N,»(v)[x'z;]  -  NO*  [x'l9*]t  0[‘s]  (66) 
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in  the  energy  region  0.1  <  KEcm  <  10.0  eV.  For  energies  below  ■'<  0.1  eV , 
the  rate  coefficient  for  the  reaction  shown  in  Eq .  (66)  is  very  small 
('^  10-1/  cm-5/ sec)  compared  to  the  I.angevin  rate  and  exhibits  a  negative 
energy  dependence  (Refs.  91-95).  A  mechanism  for  this  low  energy  behavior 


was  first  discussed  by  Ferguson  (Ref.  93)  and,  more  recently,  by  Pipano 
and  Kaufman  (Ref.  80)  and  Hopper  (Ref.  82). 


Using  our  minimum  energy  ^A"  potential  surface  (>  =  4_>°) ,  calculations 
of  the  cross-section  for  reaction  (66)  were  carried  out  employing  the  R- 
matrix  propagator  technique  (Ref.  62).  For  each  value  chosen  for  the 
collision  energy,  a  sum  over  all  allowed  j  states  was  carried  out  to 
convergence.  For  the  range  of  collision  energies  considered  here  (0.1 
<  KEcm  <  10.0  eV) ,  this  sum  converged  for  the  off-diagonal  elements  of 
the  S-matrix  for  j  <  400.  The  effect  of  increased  vibrational  energy  in 
the  N£  molecule  was  studied  by  analyzing  the  collisions  occurring  higher 
up  along  the  reaction  surface  in  the  0+  +  N2  entrance  channel. 


Our  calculated  reaction  cross-sections  are  given  in  Table  8  and 
illustrated  in  Fig.  15.  We  see  that  increasing  either  the  translational 
energy  of  the  collision  or  the  vibrational  energy  of  the  N2  molecule  results 
in  a  dramatic  increase  in  the  reaction  cross-section  at  low  energies.  For 
comparison,  the  Uangevin  rate  (Refs.  96  and  97),  which  is  determined  by 
long-range  orbiting  collisions,  is  shown  for  comparison.  We  see  that  the 
Langevin  rate  is  totally  inaccurate  both  in  magnitude  and  energy  dependence 
in  this  energy  range.  The  physical  effect  that  gives  rise  to  our  calculated 
energy  dependence  is  the  presence  of  the  barrier  in  the  reaction  channel 
which  cuts  off  the  cross  section  for  E  <  0.1  eV.  The  experimental  cross 
section  has  been  measured  by  Kaneko  and  Kobayashi  (Ref.  95)  for  E  <  0.5  eV 

where  it  eventually  turns  around  and  begins  to  increase  as  the  energy  is 
further  decreased,  reflecting  the  curve-crossing  mechanism  (Ref.  93) 
governing  the  very  low  energy  behavior  of  this  reaction.  The  experimental 
beam  data  of  Giese  (Ref.  91)  are  also  shown  in  Fig.  15.  The  overall 
agreement  with  our  ground  vibrational  state  curve  is  excellent. 


A  more  direct  comparison  can  be  made  for  the  energy  dependence  of  the 
reaction  rate  constant.  Using  an  averaged  velocity  given  by 


V, 


re  I 


2  E 


Cm 


1/2 


(67) 


the  reaction  rate  constants  have  been  calculated  as  a  function  of  the 
vibration  levels  of  N^..  These  calculated  rate  constants  are  given  in 

Table  9  and  are  illustrated  in  Fig.  16.  The  experimental  data  of  the  NOAA 

group  and  others  (Refs.  98  and  99)  are  shown  for  comparison.  The  discrepancy 

for  E  <  0.2  eV  is  again  due  to  the  fact  that  we  have  not  included  in  our 

analysis  the  dominant  curve-crossing  mechanism  for  these  low  energy  collisions. 


A  direct  comparison  of  the  relative  roles  of  vibrational  and 
translational  energy  is  shown  in  Fig.  17  where  the  vibrational  energy 
has  been  incorporated  into  an  equivalent  total  energy  for  the  reactant 
channel.  It  is  apparent  that  no  significant  differences  arise  from  the 
distribution  of  the  total  energy  into  a  vibrational  rather  than  a 
translational  mode.  The  data  of  Schmeltekopf  et  al.  (Ref.  100)  are  shown 
for  comparison,  where  the  relationship  E  ^  =  3/2(k'i'v^j))  has  been  used. 
Again,  the  discrepancy  below  'u  0.2  eV  arises  from  the  absence  of  the 
dominant  low  energy  curve-crossing  mechanism  from  our  analysis. 

A  semi-empirical  model  for  the  vibrational  dependence  of  the  0+  +  Nj 
reaction  has  been  proposed  by  O'Malley  (Ref.  101).  This  model  predicts 
that  vibrational  excitation  of  the  N2  molecule  should  affect  the  rate 
constant  much  more  strongly  than  the  kinetic  energy  of  the  collision. 

This  is  in  disagreement  with  both  experiment  and  our  calculations,  but 
the  overall  qualitative  behavior  (increase  with  v)  of  the  rate  constant 
is  exhibited  by  O'Malley's  model. 

Finally,  we  examine  the  reaction: 

0*  [4S]  +  N2[a3S^]-N20+-  products  (68) 


As  indicated  in  Fig.  2,  this  reaction  can  begin  on  the  surfaces  ^A” ,  ^A"  and 
A"  in  the  general  bent  geometry  conformation.  However,  these  excited  hyper¬ 
surfaces  surely  exhibit  crossing  with  surfaces  of  the  same  symmetry  that 
lie  at  lower  energies.  Thus,  the  reaction  given  by  Eq.  (68)  may  start  on 
the  adiabatic  surface  shown  in  Fig.  2  but  the  reaction  diabatically  can 
connect  to  several  different  product  channels.  This  situation  is  in 
contrast  to  the  ground  state  reaction  which  has  only  one  accessible  potential 
energy  surface  for  low  energy  collisions.  The  reaction  given  by  Eq.  (68) 
should  therefore  proceed  at  near  Langevin  rate  unless  an  activation  barrier 
is  present  early  in  the  reactant  channel.  Detailed  rh_  initio  studies  of 
this  hypersurface  would  appear  to  be  very  complex  and  time  consuming  since 
the  adiabatic  channel,  for  example,  is  represented  by  the  eighth  excited 
state  for  ^A"  symmetry. 


In  summary,  the  0+  +  N2  *  N0+  +  N  reaction  has  been  shown  to  take  place 
on  the  lowest  ^A  hypersurface  for  this  system,  for  collisional  energies  greater 
than  about  0.1  eV.  This  confirms  the  earlier  results  reported  by  Hopper  that 
a  bent  configuration  of  the  nuclei  exhibits  the  lowest  energy  for  this  system. 

A  small  barrier,  located  symmetrically  between  the  entrance  and  exit  channels, 
gives  rise  to  the  experimentally  observed  energy  dependence  for  this  reaction 
for  KKcm  >  0.1  eV.  The  energy  dependence  is  also  shown  to  be  nonspecific, 
either  vibrational  or  translational  energy  is  effective  in  promoting  the 
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reaction.  Finally  the  very  low  energy  behavior  (K.EC  < . 1  eV)  o i  this  reaction  must 
be  controlled  by  a  separate  reaction  mechanism  involving  dypersurf ace  crossings 
and  weak  spin-orbit  coupling  as  first  proposed  by  Ferguson  (Ref.  93).  For  such 
low  collisional  energies,  the  barrier  in  the  ^A"  surface  prevents  direct  charge 
rearrangement  collision  from  occurring  and  an  indirect  mechanism,  based  on 
doublet-quartet  surface  crossings,  must  come  into  play. 


t 
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Fig.  2.  Molecular  Correlation  Diagram  for  the  Quartet  States  of  NjO*  in  Ct  Symmetry 
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Table  2 
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Long-Range  Force  Interaction  Parameters 
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Table  4 


Long-Range  Interaction  Potentials,  0+  +  (hartrees) 


V  (x,  8,  R^  =  2.0742  bohrs) 


x  (bohrs) 

3=0° 

6=45° 

8=90° 

20.0 

-.109(-3) 

-• 553 (—4) 

-. 151 (-5) 

15.0 

— . 299 ( —3) 

-. 164  (-3) 

-. 284 (-4) 

10.0 

-.129 (-2) 

-. 781(-3) 

-.268 (-3) 

8.0 

-.296  (-2) 

- . 187 (-2) 

-. 788(-3) 

7.5 

-. 377 (-2) 

-. 242(-2) 

- . 107 (-2) 

7.0 

-.490(-2) 

-. 318 (—2) 

-.147 (-2) 

6.5 

- . 650 (-2) 

-.428(-2) 

-. 206(-2) 

6.0 

-.884 (-2) 

-.591(-2) 

-. 298 (-2) 

5.5 

-.124(-1) 

- . 842(-2) 

-.442 (-2) 

5.0 

-.180 (-1) 

-.124(-1) 

-.681(-2) 

4.5 

-. 274 (-1) 

-. 192(-1) 

-.110(-1) 

4.0 

-.440(-l) 

-. 314(-1) 

- . 187 (-1) 

5.4868 

(6=0°) 

-.1250(-1) 

[-.1508(-l)]a 

5.4740 

(6=8.0975°) 

-. 1245(-1) 

[— - 1475(— 1) ] 

5.3722 

(8=24.4103°) 

-. 1210 (-1) 

[-.1217 (-1) ] 

Ref.  78.  These  bracketed  values,  which  are  derived  from  an  MC-SCF 
calculation,  appear  to  overestimate  the  polarization  well  for  the 
linear  interaction  and  exhibit  an  incorrect  angular  dependence  for 
the  long-range  interaction. 


Table  5 


Long-Range 

Interaction 

+ 

Potential,  N  +  NO 

(hartrees) 

V  (x,  0,  RNq+  =  2.0261  bohrs) 

x  (bohrs) 

8=0° 

8=45° 

8=90° 

20.0 

-. 236(-4) 

-. 235(-4) 

-. 235 (-4) 

15.0 

-•  757 (-4) 

-.752 (-4) 

-. 748(-4) 

10.0 

-. 398 (-3) 

-.393(-3) 

-.388(-3) 

8.0 

-.101(-2) 

-. 990(-3) 

— . 969 (—3) 

7.5 

- . 133 (-2) 

-. 130(-2) 

-. 127 (-2) 

7.0 

-.178(-2) 

-. 173 (—2) 

-.169(-2) 

6.5 

-.244 (-2) 

-. 237 (-2) 

-.230(-2) 

6.0 

-. 345(-2) 

-. 333(-2) 

- . 322 (-2) 

5.5 

-. 505(-2) 

-.485 (-2) 

-.465(-2) 

5.0 

-. 769(-2) 

-.734(-2) 

-. 700(-2) 

4.5 

-.123(-1) 

-.117 (-1) 

-.110(-1) 

4.0 

-. 212(-1) 

-. 199(-1) 

-.185(-1) 

7.3298  (6=0°)  -.146(-2)  [-. 832(-3) ]a 


Ref.  78.  The  bracketed  value  reported  from  an  MC-SCF  calculation 
underestimates  the  polarization  effects  in  the  exit  channel. 
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Table  6 

Potential  Energy  Surfaces  for  N20+ 
in  the  Vicinity  of  the  Saddle  Region 

0+  +  N2 


x  (bohrs) 

V(x) 

(hartrees) 

V 

in 

CO 

0. 

+.0713 

4,5 

— 

+.0899 

4.3 

— 

+.1140 

4.2 

+.0207 

— 

4.1 

+  .0265 

j  +.1394f 

/  +.1505 ) 

4.0 

+.0633 

_ 

3.9 

+.0737 

+  .1198 

3.8 

+.1281 

_ 

3.7 

+.1475 

+.0992 

3.6 

+.2188 

— 

3.5 

— 

+.0914 

3*3 

— 

+.1166 

3.2 

— 

+.1502 

N  +  NO+ 


00 

-.0401 

+  .1536 

4.3 

4.2 

+.0530 

+.0586 

4.1 

4.0 

+.0901 

+.0603 

3.9 

3.8 

+.1278 

+.0669 

3.7 

3.6 

+.1623 

+.0852 

3.5 

3.4 

:: 

+.1213 

3.3 

— 

+.1879 
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Table  7 


A  ^ 

Reaction  Path  Potential  Energy 
for  0+  +  N2  -*•  N0+  +  N 

Collision  Angle  =  45° 


0+  +  n2 

x  (bohrs) 

00 

4.1 

3.9 

3.7 

3.5 

3.3 

N0+  +  N 

x  (bohrs) 

00 

4.1 

3.9 

3.7 

3.5 

3.3 


V  (x)  (hartrees) 

0. 

-.0050 

-.0029 

+.0037 

+.0129 

+.0313 

V  (x)  (hartrees) 

-.0401 

-.0179 

-.0102 

+.0024 

+.0240 

+.0486 


It 


Table  8 


Reaction  Cross 

KECm  CeV) 

Section  For 

0+  +  N2(v)  -*■  NO+  + 

a  (P) 

N 

v  *=  0 

v  =  1 

v  =  2 

v  =  3 

v  =  4 

.0816 

.00069 

.102 

.499 

1.20 

2.27 

.163 

.00937 

.178 

.755 

1.53 

2.53 

.272 

.0576 

.401 

1.01 

1.97 

2.87 

.544 

.339 

.876 

1.61 

2.28 

3.21 

.816 

.704 

1.27 

1.92 

2.93 

3.42 

1.63 

1.90 

2.57 

2.92 

3.74 

4.01 

2.72 

3.19 

3.46 

3.90 

4.47 

4.63 

5.44 

4.28 

4.64 

4.88 

5.01 

5.13 

8.16 

4.70 

4.81 

4.94 

5.09 

5.25 
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Table  V 


Reaction 

Rate  Constants 

For  0+  +  N2^(v) 

->  N0+  +  N 

KEcm  (eV) 

k  (cmVsec) 

v  =*  0 

v  =  1 

v  -  2 

v  *  3 

v  =  4 

.0816 

9.09(-15) 

1.27 (-12) 

6.20C-12) 

1.49(-11) 

2.82(-ll) 

.163 

1.65(-13) 

3.13(-12) 

1.33(-11) 

2.70(-ll) 

4.28(-ll) 

.272 

1.31(-12) 

9.10(-12) 

2.30(-ll) 

4 .47 (-11) 

6.5K-11) 

.544 

1.09(-11) 

2.81(-11) 

5.19(-11) 

7.31(-11) 

1.03(-10) 

.816 

2 .77 (-11) 

5.00(-ll) 

7.54(-ll) 

1.15(-10) 

1.34(-10) 

1.63 

1.06(-10) 

1.43(-10) 

1.63(-10) 

2.08(-10) 

2.23(-10) 

2.72 

2.29(-l0) 

2.48(-10) 

2.80(-10) 

3.21(-10) 

3.33(-lO) 

5.44 

4.35(-l0) 

4.71(-10) 

4 ,95(-10) 

5.09(-10) 

5.2K-10) 

8.16 

5.84(-10) 

5. 99 (-10) 

6.15(-10) 

6.33(-10) 

6.53(-10) 
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